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Abstract: 

Studies  of  the  dynamic  stabi lity  of  heavy  bodies  with  lift¬ 
ing  surfaces  towed  behind  an  airplane  are  presented.  Following  a 
fairly  complete  analysis  of  cable  dynamics,  analyses  of  the  longitu¬ 
dinal,  lateral  and  coupled  motions  of  the  system  are  developed  and 
discussed.  A  current  configuration  is  chosen  as  a  starting  point 
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NOTES:  ! )  "Upwind  wave"  and  "downwind  disturbance"  both  denote  a  wave 

traveling  upwind;  similarly,  "downwind  wave"  and  "upwind 
disturbance"  both  denote  a  wave  traveling  downwind. 

2)  Appendix  B  "T  "  in  t  is  section  is  not  as  mentioned  earlier, 
but  is  rather  Tf e  a-tual  cable  tension  at  the  vehicle  mount. 
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SUMMARY 

Heavy  bodies  with  lifting  surfaces  are  towed  behind  airplanes 
in  order  to  isolate  certain  types  of  instrumentation,  such  as  magneto¬ 
meters,  from  the  airplane.  The  imoortance  of  the  stability  of  these  veh¬ 
icles  has  recently  been  magnified  due  to  great  increase  in  the  sensitivity 
of  such  instrumentation. 

The  coupled  motions  of  the  complete  system,  which  consists  of  the 
towed  vehicle,  the  cable,  and  the  stable  tow  airplane,  are  very  complex. 
This  complexity  demands  the  use  of  certain  simplifications  in  the  analyses 
which  also  aid  the  understanding  of  the  solutions  as  well.  Different  sim¬ 
plifying  assumptions  were  studied  in  order  to  determine  those  applicable 
to  the  attainment  of  useful  results. 

A  fairly  complete  analysis  of  cable  dynamics  is  included  in  the 
study.  This  analysis  leads  to  the  primary  simplification  that  the  cable 
moves  from  one  equi librium  shape  tc  the  next  whi le  the  towed  vehicle  moves 
dynamically.  A  towed  vehicle  which  is  believed  to  have  good  characteris¬ 
tics  is  chosen  as  a  starting  point.  The  longitudinal  and  lateral  direc¬ 
tional  dynamics  are  then  studied  both  separately  and  coupled,  and  the  im¬ 
portant  parameters  affecting  the  stability  of  the  vehicle  are  isolated. 

The  results  are  presented  in  the  form  of  stability  boundaries  of  the 
dominant  parameters. 
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NOMENCLATURE 
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velocity  of  a  wave  transmitted  along  a  cable  in  a  vacuum 
drag  coefficient  of  rope 

length  of  cable  whose  drag  =  T  when  normal  to  the  flow 

center  of  gravity 

drag  of  body 

diameter  of  rope 

exponent i a  I 

aero  force/ length  to  inclined  cable 
external  force  vector 

gravitational  constant  32.2  ft/sec/ sec 

moment  of  momentum 

moment  of  inertia  about  the  x-axis 

product  of  inertia  about  y-axis 

moment  of  inertia  about  the  y-axis 

moment  of  inertia  about  the  z-axis 

imaginary  vector  notation  ~  ft 

kinetic  energy 

wave  number 

external  moment  about  the  x-axis  cable  length 

projection  of  cable  in  lateral  plane  cable  length 

external  applied  moment  vector 

external  moment  about  the  y-axis 

mean  aerodynamic  chord 

mass 

external  moment  about  the  z-axis 
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i  i  i 

projection  of  cable  in  vertical  plane 

roll  velocity  about  x-axis 

potential  energy 

perturbation  roll  velocity 

pitch  velocity  about  y-axis 

perturbation  pitch  velocity 

drag  per  unit  length  of  cable 
yaw  velocity  about  z-axis 

perturbation  yaw  velocity 

distance  measured  along  cable 
wing  area 

Laplace  transform  operator 
cable  tension 

cable  tension  rolling  moment  contribution 
cable  tension  pitching  moment  contribution 
cable  tension  yawing  moment  contribution 
cable  tension  component  along  the  x-axis 
cable  tension  component  along  the  y-axis 
cable  tension  component  along  the  z-axis 
velocity  along  the  x-axis 
perturbation  velocity  along  the  x-axis 
veloci ty 

velocity  along  y-axis 

total  velocity  of  towed  body 

freestream  velocity 

perturbation  velocity  along  y-axis 

weight  of  body 
velocity  along  the  z-axis 


wei ght/ length  of  cable 

perturbation  velocity  along  the  z-axi s 

space  fixed  inertial  axes 

body  fixed  axes 


cable  weight  to  drag  ratio  parameter  ~  tan-1  (W/2R) 
vehicle  angle  of  attack 

vehicle  sideslip  angle 

wavelength  of  traveling  wave  along  cable 
f I i ght  path  ang  le 

phase  angle 

frequency  of  traveling  wave  along  cable 
non-dimensional  height  of  a  point  on  the  cable  ~  y/c 

angle  between  cable  and  relative  wind  (See  Section  U) 

angle  between  the  cable  and  the  relative  wind 

perturbation  angle  between  cable  and  relative  wind  (See  Section 

ED 


angle  between  cable  and  body  centerline 
mass  density  of  tow  cable 

non-dimensional  horizontal  distance  of  a  point  on  the  cable  from 
the  origin  ~  x/c 

air  density 

non-dimensional  distance  along  the  cable  from  the  origin  ~  s/c 

period  of  traveling  wave  along  cable 
cable  tension  ratio  T/T0 

angle  of  inclination  of  the  cable  from  the  local  vertical 
vehicle  rol I  angle 

vehicle  yaw  angle 

general  symooi  for  frequency  ~  radians/second 


Subscri pt 


equilibrium  condition 
zero  I i ft 

Other  Notation 

second  derivative  (total  d2/dt2  or  partial  S3/St3  with  respect 
to  time 

first  derivative  (total  d/dt  or  partial  3/St)  with  respect  to 
time 

first  partial  derivative  with  respect  to  a  space  variable 
second  partial  derivative  with  respect  to  a  space  variable 


I  INTRODUCTION 


The  possibi lity  of  utilizing  airborne  towed  vehicles  as  a  solu¬ 
tion  to  a  number  of  problems  has  existed  for  many  years.  These  vehicles 
have  ranged  from  simple  weights  or  drogues  at  the  ends  of  antennas  or  re¬ 
fueling  I  i nes,  through  more  complex  machines  containing  measurement  instru¬ 
mentation  that  it  is  desired  to  isolate  from  the  main  aircraft,  to  large 
fuily-manned  gliders.  The  present  study  is  limited  to  unpiloted  systems 
and  is  primarily  oriented  toward  consideration  of  the  problem  of  the  sta¬ 
bility  character i sti cs  of  vehicles  uti  lizing  aerodynamic  surfaces. 

There  have  been  a  number  of  investigations  of  particular  aspects 
of  the  towed  vehicle  stability  problem.  For  the  most  part,  these  investi¬ 
gations  have  been  concerned  with  either  very  light  towed  bodies,  such  as 
antennas,  where  the  problem  is  primari ly  one  of  cable  dynamics;  or  with 
heavy,  high  drag  bodies  whose  inertia  prevented  the  development  of  sta- 
bility  problems.  Seldom,  however,  has  the  problem  of  a  towed  vehicle  with 
lifting  surfaces,  which  is  light  enough  to  have  appreciable  dynamic  response, 
but  heavy  enough  to  have  a  predominant  effect  on  the  cable  dynamics,  been 
considered.  In  addition,  the  increasing  complexity  and  sensitivity  of  the 
instruments  carried  in  the  towed  vehicles  has  made  the  attainment  of  some 
insight  i n+o  the  important  design  parameters  of  the  vehicle  imperative. 

The  study  reported  here  is  an  attempt  to  correlate  the  aspects  of  both 
cable  and  vehicle  dynamics  into  a  complete  system  analysis  designed  to 
isolate  the  important  parameters. 

The  total  system  consists  of  the  vehicle  dynamics,  the  cable 
dynamics,  and,  of  course,  those  of  the  mother  airplane.  It  has  been  as¬ 
sumed  from  the  outset  that  the  dynamic  properties  of  the  mother  airplane 


are  well  known  and  that  it  is  stable.  Consequently,  the  system  of  primary 
interest  is  that  of  the  combination  of  cable  and  vehicle  dynamics.  This 
system  is  depicted  in  block  diagram  form  in  Figure  I. 

The  simplicity  of  the  figure  belies  the  actual  complexity  of 
the  system.  The  input  vA  in  Figure  I  is  the  input  to  the  cable  from  a 
disturbance  of  the  mother  airplane.  This,  then,  incites  a  cable  vibration, 
m,  which  incites  the  towed  vehicle.  In  addition,  the  towed  vehicle  is 
also  subjected  to  extraneous  gust  loads  un  which  may  well  be  the  primary 
excitation.  The  combined  response  of  the  vehicle,  Vr,  which  is  a  set  of 
a  number  of  response  variables,  is  then  coupled  with  the  original  cable 
excitation  to  form  the  summed  input  to  the  cable  dynamics  is. 

The  complexity  of  the  system  would  make  a  direct  attack  on  a  com¬ 
plete  system  analysis  so  involved  that  the  possibility  of  concealing  some 
of  the  important  parameters  would  exist.  Consequently,  the  problem  was 
approached  through  a  series  of  studies  which  started  with  very  simplified 
approximations  and  proceeded  stepwise  to  the  more  complex  cases.  Not  only 
did  this  approach  insure  that  the  more  important  parameters  were  identified 
and  investigated,  but  it  also  resulted  in  an  insight  as  to  the  amount  of 
simplification  that  could  be  applied  and  still  give  reasonable  results. 

This  report  follows  the  same  line  of  attack  to  some  degree. 
Several  studies  are  presented  which  may  appear  elementary  or  even  inap¬ 
plicable.  However,  it  is  felt  that  their  inclusion  will  lead  to  a  better 
understanding  of  the  system  and  of  the  simplifying  assumptions  utilized  in 
the  attainment  of  the  majority  of  the  results. 

A  configuration  which  had  evolved  through  years  of  experimenta¬ 
tion  was  chosen  as  a  point  from  which  the  parametric  study  could  emanate. 
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As  might  be  expected.  The  basic  configuration  displayed  very  good  charac¬ 
teristics.  The  primary  results  are  presented  in  the  form  of  stability 
boundaries  which  developed  as  the  most  convenient  method  of  expressing 
the  effect  of  variations  in  the  important  parameters.  Other  results,  such 
as  the  effects  of  tow  line  variations,  are  discussed  and  comparison  plots 
are  presented.  The  agreement  of  the  responses  calculated  for  the  basic 
configuration  with  those  observed  in  flight  indicates  that  a  degree  of 
confidence  may  be  placed  in  the  results  presented,  at  least  so  far  as  an 
indication  of  the  major  trends  to  be  expected. 


-4- 


H  CABLE  DYNAMICS 

The  dynamics  of  cables  or  strings  in  vacuums  are  presented  with 
varying  degrees  of  completeness  in  many  mathematical  texts  and  other  trea¬ 
tises  such  as  References  I,  2,  and  3.  Seldom,  however,  are  these  discus¬ 
sions  carried  to  the  completeness  required  in  a  study  of  the  dynamics  of 
a  towed  vehicle  traveling  at  speeds  of  such  magnitude  that  the  aerodynamic 
or  hydrodynamic  loads  on  the  cable  are  not  negligible.  This  section  is 
an  attempt  to  correlate  the  analyses  of  somewhat  simplified  vibrating  string 
problems  with  the  more  complex  towed  vehicle  problem. 

The  classical  vibrating  string  problem  is  given  consideration 
first  in  order  to  establish  the  basic  concepts  utilized  in  the  more  complex 
problems.  The  shapesof  a  cable  in  equilibrium  in  an  airstream  are  next 
developed  for  both  a  constant  tension  and  a  variable  tension  cable.  These 
shapes  are,  in  effect,  the  shape  of  the  fundamental  mode  of  a  cable  in  a 
vacuum  with  a  distributed  forcing  function  equivalent  to  the  aerodynamic 
force. 

These  developments  are  followed  by  the  consideration  of  a  straight 
cable  at  an  angle  tc  the  airstream  and  the  great  increase  in  the  complexity 
of  the  problem  resulting  from  the  addition  of  the  aerodynamic  forces.  This 
leads,  then,  to  the  formulation  of  the  complete  problem  and  the  applicable 
simplifying  assumptions  which  allow  the  attainment  of  useful  results. 

Dynamics  of  a  Constant  Tension.  Wei Qht less  Cable  in  a  Vacuum 

The  transverse  equation  of  motion  of  a  weightless  cable  in  a 
vacuum  may  be  obtained  from  a  consideration  of  an  element  of  a  cable  dis¬ 
turbed  from  rest  in  a  vacuum  as  shown  in  Figure  2.  The  mass  density  of 


the  string  is  p,  so  that  the  mass  of  an  element  dS  of  1  he  string  is  udS. 
If  the  string  at  rest  lies  along  tr  x-axi  of  an  x-z  orthogonal  axis 
system,  the  equation  expressi  ng  the  motion  of  the  string  i  r>  the  z-direc- 
tion  is 

p.dS  32z/&t2  =  (T  +  dT)  sir,  (X  +  dX)  -  T  sin  X 
where  T  is  the  tension  in  t:.o  firing,  and  dS  is  the  elemental  length. 
Assuming  the  line  stretch  during  the  vibration  to  be  negligible,  dS  is 
approximately  a  constant,  and  dT  =  0.  The  equation  of  motion  is  Then 

udS  32z/at2  =  T  cos  X  dX 

Since  The  line  stretch  is  small,  the  angle  X  m^st  be  small  and 

COS  X  =  _ ! _  R;  I 

y/l  +  Tan2X 

when  the  square  of  the  tangent  of  a  small  angle  is  neglected  compared  to 
I.  Also, 

dx  =  ax/ as  ds 

\  s:  sin  X  Tan  X  ^  Tan  X  =  bz 
,1  4  Tan2\ 

V 


and  si  nee  dS  Ra  dx 

dx  =  a(3z/ax)  cs  =  a5z  cs 
ax  ax2 

The  equation  of  motion  is  then 

ifz  -  i  afz 

at?  u  axc 

or,  defining  T/g.  =  a2 


n- 


Equation  H-1  is  seen  to  be  of  the  form  f  the  well  known  wave  equation. 
The  most  general  solution  was  first  developed  by  D’Alembert  and  is  known 
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to  be 


z  =  f(x  -  at)  -1-  y(x  +  at)  *  XT -2 

This  general  solution  may  be  modified  to  apply  to  any  particular  case  pro¬ 
vided  the  restrictions  due  to  the  assumptions  used  in  developing  Equation 
H-l  are  not  violated. 


Infinite  String  -  Traveling  Waves 

Consider  an  infinite  string  which  extends  from  -  ®  to  +  ®  along 
the  x-axis  with  the  conditions  that  z  =  $(x)  and  dz/dt  =  Y(x)  when  t  =  0. 
Application  of  the  initial  conditions  to  Equation  TT-2  yields: 

z (X;  0)  =  gCx)  +  f  (x)  =  $(x)  H-3a 

=  a  [g»(x)  -  f»(x)]  =  Y(x)  II -3b 


3t 


t=0 


Integrating  Equation  II -3b  gives 


g(x)  -  f(x)  =  I/a  f0  Y(£)  df  +  k 


n-4 


Combining  Equations  TT-3a  and  TT-4  yields 

x 

g(x)  =  1/2  $(x)  +  I /2a  £  Y(§)  d£  +  1/2  k  n-5a 

x 

f(x)  =  1/2  §(x)  -  I /2a  Jo  Y(§)  d§  -  1/2  k  II-5b 


Substitution  of  Equations  H~5a  and  b  into  Equation  TT-2  then  yields  the 

solution  corresponding  to  the  given  initial  conditions: 

x+at 

z  =  1/2  [fCx  -  at)  +  *<x  +  at)  +  l/a  J  Y(§)  d§]  II-6 

x-at 


The  most  useful  wave  form  in  the  study  of  cable  dynamics  is  the 
harmonic  wave.  If  the  wave  profi le  at  t  =0  is  given  by 

§(x)  =  b  cos  mx 

Then  at  any  time  t,  the  displacement  z  is 


rr-7 


z  =  b  cos  m(x  -  at) 
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The  following  parameters  of  this  wave  may,  then,  be  defined 

wavelength  y  =  2tt/ m 

period  t  =  y/a 

frequency  Tj  =  I/t 

propagation  velocity  a  =  T)v 

wave  number  k  =  \/y 

Incorporating  these  definitions.  Equation  H~7  may  be  written  as 

z  =  b  cos  2rr(kx  -  Tit) 


U-8a 


or 


2tt j  ( T|t  -  kx) 


z  =  A  e  *  n~8b 

where  A  =  b  e  if  the  wave  is  out  of  phase  with  some  reference  by  the 
angle  e,  or  A  =  b  if  there  is  no  phase  difference. 

The  kinetic  energy  of  a  vibrating  string  may  be  determined  from 
the  velocity  dz/dt  of  any  point  of  the  string.  It  is 

1%  =  1/2  J*  10,(2 ) 3 dS  =  1/2  J  nCz)2dx  H-9 

where  z  =  dz/dt.  The  potential  energy  stored  in  the  string  may  be  found 
by  considering  the  stretch  in  an  element  of  the  string  from  dx  to  dS.  The 
amount  of  work  done  is  then  T(dS  -  dx)  and  the  potential  energy  is 

PE  =  J*  T(dS  -  dx) 


or,  since  dS  =  \j dx3  +  dz2 


=  JY  ['  \j  I  +  Oz/dx)2  -1  ]  dx 

Furthermore,  since  the  slope  of  the  string  is  assumed  to  always  be  small, 
the  square  root  may  be  approximated  by 


and 


y I  +(dz/dx)2  =  i  +  1/2  Oz/3x)2 


PE  1/2  J  T (dz/dx)2dx 


rt-io 
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With  a  positive  running  traveling  wave,  z  -  f  (x  -  at),  the  energies  become 
Kg  -  1/2  n  a2  J*  (f»)2  dx  =  1/2  T  /  ( f ' ) 2  dx  H-lla 

PE  =  1/2  T  J  (f ')2dx  H-l  lb 

and  they  are  equal. 


Finite  String  -  Fixed  Ends 

When  the  ends  of  the  cable  are  fixed,  the  solution  may  be  obtained 
by  requiring  that  Equation  U-6  satisfy  appropriate  boundary  conditions 
as  long  as  they  are  linear.  If  more  complicated  situations  are  anticipated, 
it  is  perhaps  more  informative  to  consider  even  the  simple  cases  by  uti¬ 
lizing  the  method  of  separation  of  variables.**  This  method  is  described 
for  this  situation  in  References  I  and  2  and  many  other  texts  and  is  there¬ 
fore  only  included  here  in  outline  form. 

Considering  the  basic  wave  equation.  Equation  XX- I,  with  initial 


conditions 


zCx,  0)  =  £(x) 

?z/3t|t=0  =  «(.,) 


n-12 


and  fixed  and  boundary  conditions 

z  (0,  t)  =  0 
zU,  t)  =  0 

the  assumption  of  a  particular  solution  of  the  form 

z(x,  t)  =  F (x)G(t ) 

and  substitution  into  Equation  IX- 1  yields: 


ET-13 


n-u 


**0ther  methods  such  as  the  Laplace  transform  method  may  also  be  used. 
In  general,  however,  complicated  initial  or  boundary  conditions  greatly 
increase  the  complexity  of  the  applications  and  they  are  not  used  here 
for  that  reason. 
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=  j-iku  n-15 

F(x)  a2  G(t) 

The  left  side  of  Equation  H-15  is  now  independent  of  t  and  the  right  side 
is  independent  of  x.  Consequently,  they  must  both  be  equal  to  a  constant, 
say  -  p2.**  Two  equations  are  then  obtained. 

F"(x)  =  -  p2F(x)  H-I6a 

&(t)  =  -  p2a2G(t)  H-I6b 

Equations  U-I6a  and  b  are  linear  differential  equations  with 
constant  coefficients.  Application  of  Equations  EL-12  and  11-13  to  Equa¬ 
tion  EL-14  yields  the  boundary  conditions  for  Equation  H-I6a  as: 

F(0)  =  F  (A)  =  0  n-17 

The  solution  to  Equation  U-I6a  is  then 

F(x)  =  A  COS  px  +  B  Sin  px  11“ 1 8 

Application  of  the  boundary  conditions  of  Equation  H-17  then  yields: 

at  x  =  0,  F(0)  =  A  =  0  II- 1 9a 

at  X  =  JL,  FU)  =•  B  sin  pt  =  0  H-I9b 

Equation  II- 1 9b  is  satisfied  when  pA  =  nrr  where  (n  =  0,  I,  2,  ...).  Con¬ 
sequently,  the  solution  to  Equation  H-I6a  takes  the  form 

F(x)  =  B  sin  viu/A  x  11-20 

The  solution  to  Equation  H-I6b  may  then  be  obtained  as 

G(t)  =  C  cos  pat  +  0  sin  pat 
or  upon  substitution  of  p  =  n tt/X 

G(t)  =  C  cos  mra/^  t  +  D  sin  nna/A  +  11-21 

Substitution  of  Equations  H-20  and  H-21  into  Equation  TT-14  yields  solu¬ 
tions  in  two  forms. 

**J  n  general,  both  +  p3  and  0  could  also  yield  solutions.  However,  in 
this  case,  it  can  be  shown  that  these  values  of  the  constant  do  not  satisfy 
the  boundary  conditions. 


o 
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2  (x ,  t)  =  E  sin  nnx/j l  cos  rrna/i  t 
or 

z(x,  t)  =  D  sin  nnx/ l  sin  nna tl  t 

Either  of  these  forms  may  be  taken  as  the  solution  or  any  linear  conk i na¬ 
tion  may  be  taken  as  the  solution.  Therefore,  the  solutions  may  be  written: 

00 

z(x,  t)  =  E  En  sin  nnx/ £  cos  nna/ 1  t  H-22a 

n  =  I 

or 

CO 

z(x,  t)  =  E,  Dn  sin  nnx/ 1  sin  nna/ 1  t  H- 22b 

n  =  I  11 

From  equation  TT  -22a 

00 

dz  =  E  En  sin  nnx  L  “  kna  sin  kna  t  ] 
at  n=l  ill 

and  at  t  =  0 

az/at]t=0  =  0 

Therefore,  it  is  seen  from  Equation  H-12  that  solution  H-22a  would  satisfy 
the  initial  condition 

az/at|+=0  =  ycx)  =o  rr-23 

Evaluation  of  z(x,  0)  from  Equation  II-22a  yields 

z(x,  0)  =  E,  En  sin  nnx/ 1  =  $(x)  U-24 

n=l  1 

A  closed  form  solution  for  En  cannot  be  obtained  from  Equation  11-24.  How¬ 
ever,  $(x)  is  an  odd  function  of  period  11  and  may  be  expanded  in  a  Fourier 
series  of  the  form 

CO 

$(x)  =  E,  br  mx/ 1 
r= 1  1 

A  comparison  with  Equation  II-24  indicates  bp  =  En  when  r  =  n.  En  is  there¬ 
fore  obtained  by  substituting  n  for  r  in  the  half  range  formula  for  the 
coefficient  of  an  odd  function. 


I  1-25 


=  2/ i  \  $(x)  sin  nrrx/4  dx 

n 

Equation  EI-22a  is  then  the  solution  of  the  problem  with  initial  conditions 

Z  (x ,  0)  =  $(x) 

Sz/^tj-t-sO  =  0 

and  En  as  given  in  Equation  11-25. 

A  consideration  of  Equation  II-22b  reveals  that 

z  (x,  0)  =  0 

and 

dz  =  £  Dn  nrra  sin  nnx  =  Y{x)  11-26 

at  t=0  i  i 

A  similar  Fourier  analysis  of  Equation  D-26  then  yields 

.1 

nrra  Dn  =  2  J  Y(x)  sin  nnx  dx  H-27 

j e  i  0  i 

Equation  U“22b  with  0n  evaluated  from  Equation  EI-27  is  then  the  solution 
of  the  problem  with  initial  conditions 

z  (x,  0}  =  0 
dz/dt| =  Y(x) 

A  solution  of  the  problem  with  the  initial  conditions  of  Equation  11-12 
is  then  a  linear  combination  of  the  solutions  given  as  Equations  II~22a 
and  II-22b  with  the  coefficients  evaluated  from  Equations  H-25  and  TX-27. 

Forced  Vibration 

When  the  string  is  subjected  to  an  externa!  force,  the  equation 

of  motion  i s  non -homogeneous  and  has  the  form 

i  d2z  -  5az  =  F(x,  t)  n-28 

a3  St3  dx3 

where  F(x,  t)  is  the  external  force  per  unit  length  divided  by  the  tension 
in  the  string.  The  boundary  conditions  correspond  to  fixed  ends 
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zCO,  t)  =  zil,  t)  =  0  n-29 

and  the  initial  conditions  are  here  specialized  to  correspond  to  a  'tring 
at  rest 


z(x,  0)  =  §(x) 


13-30 


az/at  +=0  =  o 

for  simplicity.  The  same  procedure  would,  however,  handle  the  general  iri- 
tial  conditions  of  Equation  XI— 1 2. 

The  problem  is  no  longer  homogeneous  and  a  solution  cannot  be 
derived  from  the  sum  of  several  solutions.  Howe^i  ,  a  solution  can  be 
composed  of  a  transient  parT  plus  a  particular  i n hegra I .  The  total  solu¬ 
tion  of  Equation  H-28,  with  boundary  and  initial  conditions  as  given  in 
Equations  11-29  and  11-30,  h  is  the  form 

z(x,  r)  =  zQ(x,  t)  +  z1(xf  t;  11-31 

zQ  (x,  t)  is  the  solution  of  Equation  H-1  with  the  boundary  conditions 

z(0,  t)  »  zli,  t)  =  0 


and  initial  conditions 


z(x,  0)  =  $(x) 

8z/a+|l-0  *  0 

zj_  (x,  t)  is  the  solution  of  Equation  II-2S  with  the  boundary  conditions 

z  CO,  t)  =  zU,  t)  =0 

and  the  initial  conditions 

z(x,  0)  =  dz/5t|  =  0 

The  solution  zQ(x,  t)  was  presented  previously  as  Equation  n-22a;  only 
Zj  (x,  t),  then,  is  left  tc  be  obtained. 

The  forcing  functior.  F(x,  t)  may  be  expanded  as  a  Fourier  series 


at  any  gi ven  time. 
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F(x,  t)  =  ,E  Ak(t)  sin  k((%/a)x  11-32 

k  =  l 


Similarly,  ii  may  be  assumed  that 

00 

z,  (x,  t)  =  .  E .  Bk(+)  sin  k(ouc/a)x  11-33 

x  k=l  * 

Substitution  of  Equations  U-33  and  11-32  into  Equation  H-28  after  drop¬ 
ping  the  summation  signs  for  ease  of  notation,  yields: 

Bk(t)  +  Bk  (t)k2t%2  -  a2Ak(t)  JI-34 

with  the  initial  conditions 


Bk (0)  «  Bk (0)  =  0  11-35 

The  solution  of  the  homogeneous  portion  of  Equation  TT-34  is  then 

Bk(t)  =  Cu  cos  kui^t  +  Ck  sin  ku^t  H- 36 

1  2 

The  method  of  variation  of  parameters  then  yields  the  particular  solution 
to  Equation  H-34  which  when  combined  with  Equation  H-36  gives  the  com¬ 
plete  solution  as 


Bk(t)  =  Dk  (t)  cos  ku^t  +  Dk  (t)  sin  kc%t  11-37 

1  2 

Then,  at  t  =  0,  from  Equation  EI“35 

Bk  (0)  =  0  =  Dk  (0) 

1  rr-38 

8k(0)  =  0  =  Dk  (0)  +  k<% Dk  (0) 

Evaluation  of  the  first  and  second  time  derivatives  of  Equation  11-37,  sub¬ 
stitution  into  Equation  H -34  and  multiplication  of  both  sides  of  the  result¬ 
ing  equation  by  cos  kuj^t,  yields: 

Dki  (t)  cos2  kttfc,t  -  2ki%Dkl  (t)  sin  ku^t  cos  k<%t  +  Dkg  (t)  cos  k<%  + 
sin  k%t  +  2k(%  Dk2(t)cos2k?%t  =  o2  Ak(t)  cos  ki%t  H -39 

Substitution  of  d/dt  [D^C+Jco^kuJbt]  for  the  first  two  terms  in  Equation 
11-39  and  d/dt  [6^  (t)  sin  ko^t  cos  ku^t]  +  ku^D^  (t)  for  the  last  two 


terms;  this  becomes 
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d/dt  [Dk  {+)  cos2  +  d/dt  [6k  (t)  sin  ku^t  cos  ki%t]  +  k<%  Dk  (t)  = 

12  2 

as  Ak(t)  cos  k%t 

Integration  yields 

6k  (t)  cos2  ksuot  +  6k  (t)  sin  kuiot  cos  kt%t  +  kuot>  i\  <+)  = 

1  2  +'  3 

a2  J  Ak(T)  cos  k%  TdT 
0 

where  the  constant  of  integration  is  zero.  Division  by  cos2kt«t>t  and  sub¬ 
stitution  of  d/dt  [ok  (t)  si  n  ko^t  ]  for  thv.  resulfing  second  and  third 

2  cos  kufe  i 

terms  then  gives 


Dk  (t)  +  [  Cr2  (+)  si  n  Kaut  ]  =  as 

1  dt  cos  kujt,*  co'.2ki%t  ~h 


J  Ak(T)  COS  kuJcTdT 


Another  integration  results  in 

Dk  (t)  +  Dk  (t)  si  n  km-.  f  =  a2  J  dr1  J*  Ar(t)  cos  ku^Tdr 
1  8  co:  km.t  cosskuv>Tf  0 


t  T* 

-  I 

cos^ku^-r1 

Evaluation  of  the  right  ‘and  side  by  i ntegrat i ng  f  parts  g'>ves 

t 

Dk  (t)  +  Dk  (tj  sin  ku^t  =  a2  f  A.  (t)  s:n  ku^C'-TldT  11-40 

1  2  cos  *'%+  cor  kofct  0 

Multiplication  of  Equation  X-40  by  cos  ko^t  and  substitution  into  Equation 


11-37  gives  the  desired  com: lete  solution  of  Equation  II -34  as: 

Bk(t)  =  a2  J  Ak(T)  sir‘  ku\j(t-T)dr 


XI-41 


Substitution  of  Equation  II“4l  into  H~33  a’ter  changing  the  v-eriod  from 

2tt  to  2SL  vie  las  the  steady  state  solution  of  Equation  TJ- 28. 

00  t 

z,  (x,  t)  =  £  [a2  f  Ak  (t)  sin  akn/i  (t-T)dr]sin  k  rr/jl  x  H-42 

1  k=l  k 


From  Equation  U-32  with  the  substitution  %  =  a tt/£,  it  is  seen  that  Ak(+) 
i s  gi ven  by 

i 

Ak(t)  =  2 i&  J  F(x,  t)  sin  knfJi  x  dx 


11-43 
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This  result  holds  for  any  force  F(x,  +)  as  long  as  the  force  can  be  ex¬ 
panded  in  a  Fourier  series.  An  example  would  be  a  sinusoidal  forcing  func¬ 
tion  Fx  sin  cot  applied  at  time  t  =  0  at  some  point  x0  of  a  string  at  rest 
along  the  x-axis.  $(x)  =  ¥(x)  -  0  for  this  condition,  and  the  complete 
solution  is  given  by  Equation  11-42.  F(x,  t)  of  Equation  n~43  was  taken 
as  a  distributed  force.  Therefore,  Fa  sin  tut  must  be  considered  to  act 
over  a  small  element  Ax  of  the  string  between  x  =  x0  and  x  =  Xq  +  Ax  such 


that 

then 


Fx  s  i  n  tut  =  F0  Ax  s  i  n  tut 
F(x,  t)  =  Fx/Ax  sin  tut 


and 


Xq+Ax 

A.  (t)  =  2F1/Ax£  !  sin  tut  sin  kn/i  x  dx  H-44 

x0 


Integration  yields 

Ak(t)  =  -2£i_  sin  u>t[cos  lor  Xq(cos  kn /A  Ax-I)  -  sin  knr/jfcxo  sin  krr/i  Ax] 
AxkTT  A 

Then  the  limit  as  Ax  -*  0  gives 

Afc(t)  =  2 F1/ A  sin  kn/i  Xq  sin  tut  H-45 

Substitutions  into  Equation  H-42  then  gives  z(x,  t) . 

A  _  ^  00  A  I 

z(x,  t)  =  a  v  si  n  k t\/A  Xq  sin  kTr/£  x  J  sin  tut  sin  akn/jJ  (t—r)  dT 

A  k_l  0 

Evaluation  of  the  integral  then  yields  the  final  result 

z  (x,  t)  =  a  £  si  n  kJT  xq  sin  kTT  xftu  sin  akrr/l  t  -  aktr/jR  sin  tut]  11-46 

A  A  A  tu2  -  Cakn/£)3 

Time  Dependent  Boundary  Conditions 

Time  dependent  boundary  conditions  occur  whenever  the  ends  of 
the  string  are  not  rigidly  fixed,  but  are  restricted  in  some  manner  which 
either  varies  with  time  directly  or  is  dependent  upon  an  input  from  the 
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string.  The  situation  which  exists  at  a  variable  boundary  lies  between 
that  of  an  infinite  string  problem  and  a  fixed  boundary  problem,  for  some, 
but  not  all,  of  the  energy  is  reflected  whi le  the  remainder  is  absorbed 
by  the  boundary. 

The  method  of  separation  of  var iab  les used  previously  no  longer 
applies  when  the  boundary  conditions  are  time  dependent,  for  it  is  impos¬ 
sible  to  satisfy  the  boundary  conditions  by  adjusting  only  space  dependent 
functions.  However,  by  applying  a  procedure  developed  by  R.  D.  Mindlin  and 
L.  E.  Goodman  in  Reference  4,  it  is  possible  to  separate  the  solution  in¬ 
to  two  parts,  one  of  which  is  later  adjusted  to  satisfy  the  boundary  con¬ 
ditions  on  the  other. 

For  a  string  initially  at  rest  along  the  x-axis,  the  equation 
of  motion  in  the  absence  of  an  externally  applied  force  is  given  in  Equation 
U- 1  as 


afz  -  _l  afz  =  o  n-i 

dx3  a3  3ts 

and  the  initial  conditions  are 

z  (x,  0)  =  z  Cx,  0)  =  0 

The  boundary  conditions  considered  vary  directly  with  time  -t  each  end 
and  may  be  expressed  as 


2(0,  t)  =  Ae"JUfe+  »  gift) 
z  (£,  t)  =  Be"JUi£t  =  g3  (t) 

A  solution  of  Equation  H-l  may  then  be  chosen  as 

z(x,  t)  =  £<x,  t)  +  £  pj(t)  a j (x ) 


rr-47 


n-49 


From  a  substitution  of  Equation  11-49  into  Equation  II— I ,  it  is  seen  that 
C(x,  t)  must  satisfy  the  expression 
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2fl 

dx2 


2 

-  =  -  S  [a23j C+)  a:"(x)  -  3j (t)  a j (x) ] 

at2  i =i 


with  the  boundary  conditions 


£(0,  t)  =  0x(t) 


CU,  +)  =  3S  Ct) 


and  the  initial  conditions 

2 

C(x,  0)  =  -  2 

i  =  1 

2 

CCx,  0)  =  -  x 


“  .2,  @i(+)  at  CO) 

i  =  l 

-  .1,  3:  (+)  a;  U) 

i=l 


3,  CO)  aj  (x) 
3 1  CO)  aj  Cx) 


n-50 


n-5i 


n-52 


Consequently,  it  may  be  seen  that  if  the  a;  Cx)  are  chosen  such  that  the 
right  hand  side  of  Equations  H-51  are  zero.  Equation  Et-50  may  be  solved 
by  the  method  of  separation  of  variables.  This  requires  that  the  relations 

ax  (0)  =  1  ctj,  (0)  =  0 


n-53 


ax  U)  =  0  eta  U)  =  1 

be  satisfied.  The  aj Cx)  may  then  be  chosen  as  polynomials  in  x  or  summa¬ 
tions  of  trigonometric  or  exponential  terms  whose  coefficients  are  deter¬ 
mined  from  the  relations  11-53. 

With  the  boundary  conditions  of  Equation  U-48  certain  simpli¬ 
fications  are  attained  if  the  ajCx)  are  taken  as  summations  of  trigono- 
metr i c  terms. 


(x)  =  ax  sin  %/a  x  +  bx  cos  %/a  x 

11-54 

ctg  (x)  =  a3  sin  x  +  b3  cos  tuj/a  x 

Substitution  of  Equations  U-54  into  Equations  EI-53  yields  the  a;  Cx)  as: 


ax  Cx)  =  cos  (%/a  x  -  cot  %/a  l  sin  %/a  x 
dg  Cx)  =  esc  (w^/a  i  sin  io^/a  x 


n-55 
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Substitutions  of  .Equations  H-55  into  Equations  ET-51,  11-52  and  EI-50 
give  the  desired  boundary  conditions 

C<0,  t)  =  CU,  t)  =  0  H-56 


The  initial  conditions 

C(x ,  0)  =  -  A  cos(cut,/a  x)  +  A  cot(u^/a  A)  sinfui^/a  x)  -  B  cscCcu^/a  A) 

sinfiu^/a  x)  II-57a 

£(x,  0)  =  Ajc%[cos(a3t,/a  x)  -  co’ (ui_>/a  A)  sinfo^/a  x)J 

+  Li jcu^  cscCcug/a  A)  sinfiuj/a  x)  H~57b 

and  the  differential  equation 


a2  hlL  ~  Sfl  =  0  n-58 

3x2  St2 

The  solution  of  Equation  13-58  with  boundary  conditions  of  Equation  11-56 
and  initial  conditions  of  Equations  U~67  :na/  now  be  obtained  by  separation 
of  variables  and  is  the  combination  of  the  Deviously  presented  Equations 
U~22a  and  II -22b. 

CO 

£(x,  t)  =  I  sin  nrrx  [  En  cos  nna  t  +  Dn  sin  nna  t  ]  H-59 

n=l  l  l  i 

Substitution  of  Equation  U-57a  for  $(x)  and  Equati on  IT  -57b  for  Y(x)  in 

Equations  JIT -25  and  H-27,  and  performance  of  +Ke  necessary  integrations 

yields  the  coe-f f icien+s 


En  = 


(~n_  2B nn 


Dn=-- 


C-I)n  26iWa)_ 


n-60 


^2[(nn/i)3-  (uu^/a)2]  ’  X[(nn/X)2-  (to^/a)2] 

The  solution,  obtained  from  substitutions  of  Equations  H-60  into  Equation 
ET-59  and  the  result  along  with  the  a;  Cx)  and  |3j  (t)  into  Equation  HE-49, 
i s  then 

z(x,  t)  =  Atcos(ufc,/a  x)  -  cot(u>,/a  A)  sin(i%/a  x)]e  J(U°+  +  Be  * 


csc(coj/a  A )  sin(oj^/a  x)  + 


2  ( - 1 )  nB 


n  1  X[(nTT/X)3-  (cuj/a)2} 

-  j  u^/a  sin  (nna/i  t)]  sin(nrr/A  x) 


[rura/i  cos(nira/£  t) 


n-ei 
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The  denominator  of  the  summation  in  Equation  XT-61  is  seen  to  vanish  when 
the  end  frequency  cog  is 

u)£  =  ann/X 

If  the  corresponding  mode  (n)  is  excited,  z(x,  t)  -*  <»  and  the  frequency 
co^  is  known  as  the  n^1  modal  resonance  frequency. 

The  Development  of  the  Static  Eaui librium  Equations  of  a  Weiohtless  Cable 
Supporting  a  Body  in  a  Constant  Velocity  Flow  Field 

The  development  here  of  the  equilibrium  equations  of  a  wieght- 
less  cable  supporting  a  body  in  a  flow  field  follows  the  procedure  described 
in  Reference  5. 

The  wire  is  taken  to  be  weightless  and  i nextensi b le.  R  is  the 
drag  per  unit  length  of  the  cable  when  at  right  angles  to  the  flow. 

R  =  C0r  d  p  V3  31-62 

When  the  cable  is  inclined  at  an  angle,  X,  (see  Figure  3)  to  the  free  stream, 
the  force  per  unit  length  of  the  cable  at  right  angles  to  the  cable  is  of 
the  following  form; 

F  =  R  sin3  X  n-63 

Since  the  aerodynamic  force  on  the  cable  is  perpendicular  to 
each  element  of  length,  and  because  it  is  weightless,  the  tension  in  the 
cable  will  be  constant,  denoted  by  T,  and  has  the  following  form; 

t  scp/as  =  f  rr-64 

Where  dS  i  s  an  element  of  length,  and  <p  is  the  inclination  of  3S  from  the 
vertical  (see  Figure  3).  The  tension  may  also  be  expressed 

T  *  Rc  n-65 

where  c  is  the  length  of  wire  whose  drag,  when  normal  to  the  flow,  is  equal 
to  the  tension  T.  The  differential  equation  for  the  shape  of  the  wire 


becomes 
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dS/dcp  =  c  sec8  cp  H-66 

Next,  the  wire  may  be  considered  to  be  supporting  a  body  of  weight 

T  and  zero  drag.  The  system  is  illustrated  in  Figure  4,  where  the  origin 

0  is  taken  at  the  bottom  point  of  the  wire.  Integrating  Equation  H-66 

\ 

S  =  c  tan  cp  H-67 

The  rectangular  coordinates  of  any  point  of  the  wire  with  re¬ 
ference  to  the  origin  may  be  derived  from  the  differential  equations, 

dx/dS  =  si  n  cp  11-68 

dy/dS  =  cos  cp  ET-69 

Integrating  these  expressions  with  respect  to  the  origin  at  0 


as  before. 


S  =  c  tan  cp  =  sinh  y/c, 
x  =  c(cosh  11-1). 

Def i ni ng 


H-70 

n-7i 


cr  =  S/c,  5  =  x/c,  71  =  y/c 
Equations  11-70  and  DL-71  can  be  rewritten  as  follows 

a  =  tan  cp  =  sinh  T),  IX- 72 

5  =  cosh  7]  -I .  IT-73 

Next,  the  length  of  wire  BA  as  shown  in  Figire  5  may  be  considered.  At 
8  the  cable  tension  T  has  horizontal  and  vertical  components  capable  of 
supporting  a  body  of  weight  W  and  drag  D  provided, 

W  =  Rc  cos  cp  U-74 

D  =  Rc  si  n  cp  H-75 

The  shape  of  the  cable  supporting  a  body  may  be  determined  by 
using  Equations  11-72,  73,  74,  and  75  in  the  following  manner.  The  values 
at  point  A  are  denoted  by  subscript  I,  and  B  by  subscript  2. 
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=  sinh  %  =  tan  9  =  D/W  n>76 

a3  =  si  nh  7)g  11-77 

gx  =  cosh  %  -  I  *  sec  9  -I  H-78 

g2  =  cosh  1^3  -I  n-79 

and 

S3  "Si  -  ct  -  a/c  n-80 

ITa  -  %  =  0  =  b/c  n-81 

<7g  -  ctx  =  CT  =  s/c  n-82 


A  combination  of  these  equations  wi II  then  eliminate  the  co¬ 
ordinates  (g,  T),  a)x  and  (g,  7],  a)3  as 

ct  +  tan  9  =  sinh  (0  +  %)  =  sinh  (0  sec  9)  +  cosh  (0  tan  9)  11-83 

and 

a  +  sec  9  =  I  +  g3  ®  cosh  (0  +  U )  =  cosh  (0  sec  9)  +  sinh  (0  tan  9)  Et-84 
Thus  the  two  equations  which  determine  the  shape  of  the  weightless  cable 
are: 

a  =  sinh  (0  sec  9)  +  [cosh  (0  -  I)]  tan  9  H-85 

a  =  sinh  (0  tan  9)  +  [cosh  (0  -  I)]  sec  9  11-86 


The  development  of  the  equations  governing  the  equi librium  shape 
of  a  heavy  cable  in  an  airstream  follows  the  procedures  of  References  6 
and  7.  As  to  be  expected,  this  development  parallels  that  of  the  weight¬ 
less  cable  in  many  respects.  However,  certain  basic  differences  do  exist; 
and  since  the  resulting  expressions  describe  the  actual  physical  situation, 
they  are  developed  here. 
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The  cable  is  considered  to  have  a  weight  w  per  unit  length  and 
a  drag  R  per  unit  length  when  rerpendi cu  lar  to  the  free  stream  velocity. 

R  has  the  form, 

R  =  C0  d  p  V2  U-62 

where  C0r  is  the  drag  coefficient  of  the  cable  and  d  is  the  diameter. 

When  inclined  to  the  free  stream  at  an  angle  X  (see  Figure  6).  The  aero¬ 
dynamic  force  per  unit  length,  F,  perpendicular  to  the  cabie,  wi II  be  of 
magnitude, 

F  =  R  s>  n2  X  n-87 

For  convenience  the  development  will  be  made  first  for  a  heavv 
body  with  zero  drag  and  then  for  a  heavv  body  with  drag. 

Taking  0  as  the  origin  cf  coordinates  (see  Figure  6;  with  X  as 
defined  earlier,  the  following  relations  may  be  established, 

dx/dS  =  cos  X  11-88 

dz/dS  =  sin  X  0-69 

where  dS  is  the  element  of  length  inclined  at  an  angle  X.  The  element 
dS,  at  point  P,  is  in  equilibrium  and  thus  the  forces  along  the  element 
of  length  may  be  written 

dT/9S  =  w  sin  X  LI-90 

which,  when  integrating  S  from  the  origin  to  point  P,  gives 

T  =  T0  +  wy  II-9 I 

where  T0  is  the  weight  of  the  body  being  supported. 

Resolving  forces  perpendicular  to  the  element  of  length,  the 
second  fundamental  equation  is  obtained 


T  fo/dS  =  w  cos  X  -  R  sin3  X 


n-92 
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At  this  point  it  is  convenient  to  introduce  the  non-dimensional i - 

zing  coefficients  defined  as  follows 

T  =  T/T0,  1]  5  _z _  ,  |  =  _x _  ,  T  =  _  '  H-93 

T0/R  T0/R  T0/R 

where  T,/R  is  the  length  of  cable  whose  drag  when  perpendicular  to  the 
free  stream  is  equal  to  the  tension  T0.  Final ly,  the  weight  to  drag  ratio 
os  the  cable  p.  is  defined  as: 

w/R  =  p,  =  2  tan  a  H-94 

Uti  lizing  the  above  notation.  Equations  H-90  and  H-92  become 

t  -  I  +  p.1)  H-95 

and 

t  BX/Bct  =  u  cos  X  -  sins  \  U-96 

Elimination  of  t  from  these  two  equations  yields  the  relationship 
between  T]  and  X.  The  corresponding  values  of  §  and  j  are  found  by  evaluat- 
i  ng  the  i ntegra I s 

5  =  cot  X  311  11-97 

and 

nil 

a  =  f  cosec  X  BT)  11-98 


For  a  heavy  cable  the  process  is  as  follows;  from  Equations  11-95 


and  U-96, 


(I  +  p,T])  sin  X  B  /BT|  =  cos3  X  +  p,  cos  X  -  I  H~99 

which  in  integral  form  becomes 

nil  jr/2 

f  u3Tl  =  r  _ ti  sin  X  BX  HI- 1 00 

0  I  +p,Tj  0  i  “  p.cos  -  cos3 

From  a  substitution  of  2  tan  a  for  p,  and  evaluation  of  the  integrals,  it 


i  s  seen  that 
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log  ( I  +  2T]  tan  a)  =  sin  a  log  cos  a  +  (I  -  sin  a)  cos  X  H-101 

cos  a  -  (I  +  sin  a)  cos  X 

This  equation  was  then  used  to  calculate  T]  as  a  function  of  X.  The  coi — 
responding  values  of  §  and  a  were  obtained  by  graphical  integration  of 
Equations  ET-97  and  H-98. 

Next  the  equations  wi 1 1  be  developed  for  a  body  with  drag.  The 
towed  body  is  assumed  to  have  a  net  weight  of  W  and  a  drag  D.  Net  weight 
is  simply  the  weight  minus  any  lift  produced  by  the  body.  The  angle  that 
the  cable  makes  with  the  body  is  determined  from  the  equation 

Tan  X  =  W/D  II- I 02 

The  tension  at  thi s  point  wi I  I  be 

T  =  fFTtf  11-103 

Knowing  a  as  a  function  of  X,  the  point  of  connection  can  now  be  found; 
butT0,  the  fundamental  tension,  must  still  be  determined.  From  Equation 
EL-95,  the  value  of  T/T0  may  be  determined  at  any  point;  and  when  used  in 
conjunction  with  Equation  U-IOI,  a  graph  may  be  constructed  to  determine 
T/T0  as  a  function  of  X  directly. 

The  results  of  this  analysis  are  presented  in  graphical  form  as 
follows;  Figure  7  is  simply  the  cable  shape  in  non-dimensional  form,  the 
parameter  a  being  constant  for  each  curve;  lines  of  constant  a  are  also 
shown.  Figure  8  is  a  plot  of  X  versus  T0/T,  mentioned  earlier. 

Dynamics  of  a  Constant  Tension.  Weightless  Line  in  Moving  Air 

The  aerodynamic  force  per  unit  length  normal  to  an  element  of 
cable  at  an  angle  X  with  the  ai  r  stream  is  obtained  from  Equations  H-63 
and  U-62  as: 

dFa  =  C0  3  p  Vs  sin®  X  dx  EL- 104 
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Consideration  of  a  straight  cable  initially  situated  at  an  angle  X0  to 
the  airstream  velocity  V0  such  that  the  disturbed  angle  of  attack  is  never 
negative,  shown  in  Figure  9  where  the  x-axis  is  taken  along  the  undisturbed 
cable,  leads  to  the  following  relations  for  V  and  X. 

vs  =  cv0  -  az/at  sin  x0)2  +  caz/at)2  cos3  x0 

and  since  X0  is  fairly  large,  cos2  Xq  <  I  and  (3z/3t)2  cos2  X0  may  be  ne¬ 
glected  when  compared  with  (V0  -  3z/3t  sin  \q)3  .  Therefore 

V  *,  V0  -  az/at  sin  X0  II- 105 

From  Figure  9  it  is  seen  that 

X  —  Xq  —  X2  —  x^ 

and 


Xx  =  sin”1  az/at  cos  X.% 
V 


and  again  since  3z/3t  «  V 

Xx  az/at  cos  Xr 
V 


A  I  so, 


X2  =  tan”1  3z/3x  sin"1  3z/ax 
\3  3z/3x 


Then 


x  =  Xq  -  3z/ax  -  az/at  cos  x0/v 

or,  defining 

X  =  +  3z/3x  +  az/at  cosX0/V 
X  »  X0  -  X 

Application  of  the  trigonometric  identity 

sin2  X  =  1/2  (I  -  cos  2X)  =  1/2  (I  -  cos  2XQcos  2X  - 


U-106 


sin  2Xcsin  2X) 
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then  yields,  since  K  «  \q, 

sinaX  =  1/2  (I  -  cos  2Xq)  ~  X  sin  2X0 


or 


sinaX  =  sinaX0  -  X  sin  2X0 
or  upon  subs+i+u+ion  for  X 

sinaX  =  si naX0  “  Oz/Sx  +  9z/5t  cosX0/V)  sin  2\0  H-107 

Substitution  of  Equations  12-107  and  H-105  into  Equation  22-104,  upon 
neglecting  higher  order  small  terms,  gives  the  aerodynamic  force  as 
dFa  =  C0  p  3  V02  [sin2\0  -  sin  2X0  3z/5x  -  l/V0  (cos  X0  sin  2\0  + 

2  sin3  X0)  3z/3t]dx 

Multiplication  of  the  night  hand  side  of  this  equation  by  sin  2X0/sin  2X0 
and  introduction  of  the  constant  parameters 

5  =  Cp  p  d  Vp si  n  2X0 

T 

6  =  cos  \0  +  2  sin3  \q 

sin  2  X0 

then  gives  the  aerodynamic  contribution  in  the  form 

dFa/T  =  (§  si n2X0/si n  2X0  -  %  3z/3x  -  l/V0  ?  &  3z/3t)  dx  11-108 
Substitution  of  the  aerodynamic  force  into  Equation  22." '  yields  the  equa¬ 
tion  of  motion. 

3az  +  jr  si  n3X0  -  5  -  J_  §  <5  =  _!_  d3z 

3x2  si  n2X0  3x  V0  3t  a3  3t2 

However,  since  the  string  under  consideration  is  initially  straight  along 

the  x-axis  and  at  rest,  32t/3xa  =  3z/3x  =  3az/3ta  =  3z/3t  =  0  initially 

and  it  is  seen  that  the  constant 

§  sin3X0/sin  2 X*,  =  0** 


**"For  a  heavy  cable,  this  expression  is  f  s ? n3^o  =  £  co  s 

s  i  n2:Xo  a3 
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The  equation  of  motion  is  then 

a2z/dx2  -  5  az/ax  -  §6/v0  az/at  =  i/a2  a2z/at2  rr-109 


Forced  Infinite  Cable 

The  steady  state  transverse  response  of  a  cable  forced  by  a  sinu¬ 
soidal  forcing  function  Ae'j^  at  a  point  x  =0  will  assume  the  frequency 
of  the  forcing  function.  For  a  cable  which  extends  to  infinity  both  up¬ 
wind  and  downwind  of  the  x  =  0  point,  the  resulting  displacement  along  the 
cable  may  be  assumed  to  be  given  by  Equation  U-8b  with  T|  replaced  by  tu/2n. 

z  =  Ae  211 J  to/2n  t  -  kx)  H-IIO 

Substitution  of  this  form  into  the  equation  of  motion.  Equation  H-109, 
yields  the  desired  expression  for  k. 

k2  -  j  £/2tt  k  +  iMrF’CSauj/Vo  j  -  u/7a2)  =0  H-IM 

Then,  since  the  resulting  expression  for  k  is  complex,  it  may  be  assumed 
that 


k  =  u  +  jv 

Substitution  into  Equation  II- 1  1 1  yields 

u2  +  ?/2tt  v  -  va  -  u^^n^a3  +  j(2uv  -  §/2tt  u  +  £6u)/4tt2V0)  =  0  H-l  12 
Then  from  the  imaginary  part  of  Equation  H-II2 

2rr  v  =  S/2  (I  -  6u>/2ttV0u)  11-113 

Substitution  of  Equation  XI-II3  into  the  real  part  of  Equation  II-II2  then 
g  i  ves 

2ttu  =  +  [  I /2(  (w3/ a2  -  52/4)  +  [Cu)2/a2  -  ?2/4)2  -  68cu853/V02]l/a)  ]l/s 

H-l  14 


Then  from  Equation  H-110 

.  2nvx  j(u)t  -  2ttux) 
z  =  Ae  e 


H-l  15 
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A  downwind  disturbance  is  thus  damped  when  v  is  negative  and  amplified 
when  v  is  positive.  The  reverse  is  true  for  an  upwind  wave  which  corres¬ 
ponds  to  a  negative  value  of  u.  For  a  negative  value  of  u,  v  is  always 
positive  and  an  upwind  wave  is  therefore  always  damped.  The  damping  of 
a  downwind  wave  is  dependent  on  the  inequality 

6u)  >  I 

2ttV0u 

for  damping.  Or,  since  the  propagation  velocity  for  this  wave  is 

a  =  a)  U- 1  1 6 

2nu 

the  damping  inequality  may  be  expressed 

6  _a  >  I  H- 1  I  7 

V0 

The  term  6  is  in  the  range  f  <  6  <  00  when  Xo  is  in  its  normal  operating 
range  0  <  X0  <  90c .  Therefore,  for  very  large  X0near  90°,  the  downwind 
wave  will  be  damped  for  all  reasonable  velocities.  However,  for  Xq  near 
zero,  the  damping  depends  on  a/V0.  Then,  from  Equation  H-II4,  as 


Xq  — *  0 


and  from  Equation  U-II6 


2nu  +  \J  1/2  u)/a  Ctu/a  +  I ) 


a  -*  +  \J  2a2u)/tu+a 

Then  for  a  lightly  loaded  cable  or  high  forcing  frequency 

a  -»  +  \pl  a 

and  for  a  highly  loaded  cable  or  low  forcing  frequency 

a  \j  2auu 


Therefore,  to  insure  damping  at  low  values  of  X0  will  always  require  the 
use  of  a  light,  heavily  loaded  cable. 
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f-inite  Cable  -  fixed  Ends 

The  finite  cable  with  fixed  ends  rrav  be  obtained  from  the  pre¬ 
viously  considered  infinite  cable  dv  corsi  dering  the  cable,  at  an  angle 
X0  tc  the  ji  ret  r'ea'n,  to  be  clapped  at  two  points  along  the  x-axis,  x  =  0 
and  *  =  i.  The  equation  to  be  soivc-d  is  then  the  same  as  given  in 
Equat i on  XI- 109 

32z/3x2  -  ?  3z/3*  -  56/V0  3^/ar  =  I/a2  3sz/3t2  H-109. 

with  the  boundarv  conditions 

z(0,  t)  =  z(i,  t)  =  0  rr-i 18 

and  the  initial  conditions 


?(x,  0)  =  $(.*) 


3z/3 1 


1=0 


K>.1 


rr-i  i9 


Application  of  tfe  'netted  of  separation  of  variables  with  the 
assumed  soiutior. 

z(x,  t)  =  ?  Fk(+»  G.  (x)  11-120 

k=0  K 

where  the  G|.(rJ  are  er  t1' pgo'-.a  !  ,  leads  t0  the  t  wo  eq  eat  ions 

G"  (>  )  •  §  C»  *  naC  =  0  II-  1 21  a 

r( t)  *  a2§6/V0  F  (+ )  r  a2n2F  =0  012  lb 

A  solution  of  Fqur4  ion  IX  - 1 2 1  a  after  aop  I  i  CaJ  i  on  of  the  boundary  conditions 
of  equation  H~-i  !8  is  then 

Gk(x)  =  x  Ak  sin  kn/l  x  H-122 

and  n2  =  k2TT2/42  -  (§/2)a. 

Similarly,  a  solution  of  Tauation  U*  I2lb  is 

Fk(t)  =  e-a2§6 /2V0  +  [Bk  cos  a(-k2TT3/je2  +  ^/4(|  „  a262/4V0a)]l/3t 

h  Ck  sin  a[k2n2/ i,3  -r  ff/A  (I  -  a262/4v02)  /2t  ] 


H-123 


Substitution  of  E qudf  ion3  IJ—  122  and  i23  into  equation  XE-120,  gives  the 


resu I t 


t  C  x  ,  t )  =  c 


'■S' 


-  oa|6/2V0  t)[D^  cqs  a[k2n2/J!3  +  ?a/4(I  _  a252/4v03)  ]i/st 


+  t.  sin  a[k2rr/£s  +  |2/<i(l  -  j262/4V02)  ]1,/2+]  sin  krr/X  x 


n-124 


■<_r  ) 


J/O  . 

e*  *  A  sir-  k'Tx/i  i  €  orlhogonal,  and  is  simply 


scale  ’ac*or,  it  •'ijy  be  soecitied  fhat 


c-5-'2 


■-1  knx/jt  =  Hk  (x) 


EL-125 


vtere  H.  i  o't,_ogor>ol . 


fli  r  iicdhc.  of  tKe  ivrial  conditions  to  equation  JQ--I24  with 
rq,'t  or  JT-  i?S  incorporated  voids 

n|,  C»K  Kk(y)  =  • 

ang 

Z  [.-[k2r^/i2  +  §2/4(|  -  a262/4V02)]l/2  E.  - 
n=  I  K 

§6as/?V0  Dk]  HR  v  <  >  =  Y'x) 

T ron  vih:ch  the  coe'f  icie^'s  Dk  and  E  |,  'nay  be  determined  as 

Dk  -  <£l  dx 


p  Hk(x)2  d* 

°o 

J  o  ^ 

L  _  Jo  Y( xi  r.^C-')  dx  *  \bc?/2V0  ^  $(x)  Hk(x)  dx 


U-  126a 

n-l26b 


[x2TT2/i2  «■  §24(  I  -  a262/4V, 


Wo2)]1/2  p  Hk(x): 
J  n 


dx 


The  so !  u r !  o  -  ■  or  Equation  H-IG9  v/i  fh  the  ;ni  +  ia!  and  boundary  conditions 
of  Equations  Jj[-M8  and  '19  is  then  given  by  Equation  IX—124  with  the  coef- 
ficients  eva'uated  f^onr  Equations  Il--I26a  and  I26d„ 


Time  Dependent  Boundary  Conditions 

Consideration  of  the  same  time  dependent  boundary  problem  con 
sidered  previously  for  a  line  in  a  vacuum,  with  the  added  aerodynami c 
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parameters,  indicates  the  increase  in  complexity  created  by  the  aerodynamic 
terms.  The  equation  of  motion  is  given  in  Equation  U-109.  The  boundary 
conditions  are: 


z  (0,  t)  =  Ae  (t) 

z(X,  t)  =  Be"Ju,4+  =  g2  (t) 
and  the  initial  condi i ions  are: 

z(x,  o)  =  az/at  |+=0  =  0 

The  same  procedure  used  in  the  previous  time  dependent 
is  followed.  A  solution 

2 

z(x,  t)  =  y (x,  t)  +  E  3;  (+)  a,  (x) 

i  =  i 


n-127 


II- 128 

boundary  problem 


II- 128 


i.s  assumed,  where  y(x,  t)  must  satisfy  the  equation 


a3  32y/3x2  -  a3  ?  3y/3x  -  a3  §6/V0  3y/3t  -  a3y/at2  =  .E-a2^  (t)aj  "(x) 


+  a3  ?  3j  (t)  ai’(x)  +  a3§6/V0  $j#(t)  a;  (x)STf  (t)a;  (x) 
with  the  boundary  conditions 


n-129 


2 

y (0,  t)  =  3x(t)  -  .Ej  3; (t)  a; (0) 


n-130 


y (X,  t)  =  e8(t)  -  ,E(  8j  (+)  aj  CX) 


and  the  initial  conditions 


y  (x,  0)  =  -  E  0  i  (0>  a;  (x) 
i  =1 


3y/3t 


t=0  =  "  j?,  M°> 


n-i3i 


The  aj  (x)  are  then  chosen  to  make  the  right  hand  side  of  Equations  U-130 
van i sh. 

(x)  =  cos  u^/a  x  -  cot  i%/a  X  sin  cu^/a  x 
02  (x)  =  esc  w^/a  X  sin  tu^/a  x 


H-132 
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SubsM+u+ion  into  Equation  H-129  yields  the  equation  of  motion  in  the 
form 

a2  a3y/dx2  -  a2f  3y/dx  -  a2§6/V0  }y/3t  -  d3y/St2  = 

a5[“''1toAe  rx  (x)  +  uj£  esc  ajjj/a  l  e  r2(x)]  U-133 

where 

^  (x)  =  (cot  ui-j/a  l  +  j  a6/V0)cos  Ofc,/a  x  +  (I  -  j  a6/V0  cot  u^/a  Z)  • 

sin  cut>/a  x 

and 

r2(x)  =  cos  u j^/a  x  -  j  a6/V0  sin  ca^/a  x 
Application  of  the  method  of  separation  of  variables  leads  to  an  assumed 
solution  of  Equation  II- 1 33  of  the  form 

y(x,  t)  =  |  Xn(x)  Tn(t)  n-134 


and  Xn  is  assumed  to  be  an  orthogonal  function.  The  rx (x)  and  r8(x)  may 
then  be  expanded  in  series 

rl  =  n  =  l  ^n  ^n 


r2  (x)  -  £  Hsn  Xn 

n=l  n 


n-135 


where  H, 

1  n 


and  Hgn  are  given  by  the  orthogonality  relations 
Hin  =  ^  ri<x>  Xn  dx/  ^  Xn3  dx 

^n  =  Jq  rs(x)  Xf1  dx/  £  *n  dx 
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The  boundary  conditions,  from  Equations  11-132  and  130,  are 

y(0,  t)  =  y (JL,  t)  =  0 

and  Tne  initial  condi tions, from  Equations  EL-131  and  132,  are 


H-137 
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-(x,  0)  =  -A(cos  %/a  x  -  cot  %/a  l  sin  %/a  x)  -  B  esc  u)^/a  i  • 


sin  (u 


i/a  x 


dy/dt 


t  =0 


*  ~  n-138 

=  Aj%  (cos  (%/a  x  -  cot  (%/a  i  sin  %/a  x)  +  Bjtu^  esc  (0^/ a  l  * 


si  n  ui^/a  x 


Substitution  of  Equations  H-134  and  135  into  Equation  IE-133,  separation 
of  the  variables,  equating  the  space  and  time  dependent  portions  to  the 
constant  -  yn3  and  satisfying  the  Doundary  conditions  yields  the  space 
dependent  solution. 


Xn  =  Cn  e^x  sin  m/ Z  x 


IE-139 


the  auxiliary  equation 

yn2  =  a2[(nn /  l)a  +  (f/2)2] 
and  the  time  dependent  solution 

Tn  =  e&t  [Dn  cos  / Yn2  "  Zbf  +  En  sin  /Yn2  "  Cb  t]  +  Tp  11-140 
Tp  is  the  particular  integral  resulting  from  the  right  hand  side  of  Equa¬ 
tion  IT- 1 33  and  Cb  =  a256/2Vc.  Tp  may  be  determined  by  variation  of  para¬ 
meters  to  be 


TP  - 


_aL 


\/an2-Cb2  Yn2“two2-^Jab  Cb 


Hi  n  [  yvn2-Cb3  e'J0Jo+  -  e"Cb+(  • 


^os  An2"£b2  +  ‘  sin/  Vn3-Cb2  +  >3  -  ■  - —  - 

Yn  ~2^b 

[  /Yn2-Cb2  e"^f  -  e"£b+  (  /  YnS-Cb8  cos/y^C,2  t  -  (Cb-j»i>  ’ 

sin  /Yn2-cb2  t)  ]  H-I4I 


The  parameters  H1  n  and  Hg  may  be  obtained  from  Equations  U-136.  They 
are  complex  constants  in  yn>  £b  and  5  and  are  very  lengthy  and  are, 
therefore,  not  presented  since  they  would  only  complicate  the  picture  with 
out  adding  any  important  parameters.  They  both  contain  the  factor  l/Bn 
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and  the  quantity  Tp*  =  Bn  Tp  may  therefore  be  defined  without  loss  of 
generality.  The  solution  of  Equation  H-133  may  then  be  expressed 


V(x,  t) 


g  &(C/2  x  -  Cb1^ 


[  Fn  cos  fan  ~  ^b 
si n  nn/i  x  +  Tp*  e 


t  +  Gr 


§/2  x 


si  n 


sin  fan~^b  +3  * 

m\/l  x  TT-  1 42 


where  the  constants  Fn  and  Gn  are  obtained  by  satisfying  the  initial  condi¬ 
tions.  The  final  form  of  z(x,  t)  could  then  be  obtained  by  a  substitution 
into  Equation  H-128. 

The  example  considered  here  is  that  of  a  cable  forced  at  both 
ends.  Energy  is  supplied  to  both  ends  of  the  cable  at  constant  frequen¬ 
cies  and  the  resulting  wave  form  is  that  of  a  standing  wave.  The  waves, 
however,  are  damped  with  both  time  and  di s+ance  upstream  and  amplified  with 
distance  downstream.  If  the  standing  wave  is  considered  as  a  series  of 
direct  and  reflected  traveling  waves,  then  all  waves  introduced  at  the 
downstream  end,  whether  direct  or  reflected,  are  damped.  However,  waves 
introduced  at  the  upstream  end  may  be  damped  or  undamped  depending  on  the 
distance  of  the  point  in  question  from  the  initial  end  and  the-  wave  velo¬ 
city.  The  time  required  for  the  damping  to  override  the  afup  I  i  f  i  cat  i  on 
at  any  given  point  along  the  cable  is  then 

td  =  V0x/ a2  6 

and  the  primary  items  of  importance  are  again  the  freestream  velocity,  the 
angle  of  the  cable  with  the  stream,  and  the  tension  to  mass  ratio  of  the 
cable  as  discussed  previously. 

The  end  conditions  representing  the  case  of  a  towed  vehicle  are 
somewhat  different  than  in  the  example  considered.  The  end  conditions 
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for  the  towed  vehicle  would  actually  be  the  dynamic  equations  of  motion 
of  the  system.  Energy  is  then  both  added  and  removed  from  the  cable  dur¬ 
ing  different  periods  of  time.  The  motion  of  vehicles  of  this  type  are 
approximately  sinusoidal  and  the  case  just  considered  would  be  a  fairly 
good  approximation  of  the  energy  input  contri ubt ion.  However,  the  en¬ 
ergy  removed  from  the  string  disturbs  the  reflexion  relations  in  both  am¬ 
plitude  and  phase.  Consequently,  cancellation  of  the  downwind  wave  does 
not  occur  in  the  same  manner.  If  the  energy  inputs  and  removals  occurred 
at  the  same  frequency,  a  standing  wave  simular  to  the  example  would  exist. 
However,  they  are  not  and  the  resulting  wave  form  is  that  of  a  traveling 
wave  superposed  on  a  standing  wave. 
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IH  EQUATIONS  OF  MOTION  OF  A  TOWED  BODY 

The  equations  of  motion  of  a  towed  body  are  essentially  the  same 
as  the  dynamic  equations  expressing  the  motion  of  a  free  body  such  as  an 
airplane.  The  primary  difference,  aside  from  possible  configuration  dif¬ 
ferences,  is  the  pull  of  the  tow  cable  which  varies  with  the  dynamic  mo¬ 
tion  of  the  towed  body. 

This  section  will  briefly  describe  the  general  development  of 
the  dynamic  equations  of  the  towed  vehicle  and  develop  useful  representa¬ 
tions  of  the  cable  force. 

General  Development 

The  general  development  of  the  basic  equations  of  motion  of  a 
vehicle  in  air  are  well  described  in  References  8,  9,  and  10.  Consequent¬ 
ly  only  a  brief  outline  will  be  presented  here. 

Figure  10  demonstrates  the  configuration  considered.  The  x,  y, 
z  axis  system  is  fixed  at  the  center  of  gravity  of  the  body.  The  X,  Y, 

Z  axes  are  fixed  in  space.  Application  of  NewtonTs  second  law  then  gives 
the  vector  equations  of  motion 

F  =  m  dV/dt  DU-  I 

M  =  dh/dt  EH-2 

where  F  and  M  are  the  external  applied  forces  and  moments,  V  is  the  velo¬ 
city  of  the  center  of  gravity  of  the  body  of  mass  m,  and  h  is  the  total 
moment  of  momentum  of  the  body.  These  equations  are  referred  to  the  X, 

Y,  Z  axes  fixed  in  space.  If  these  axes  are  considered  to  be  fixed  at  the 
center  of  gravity  of  the  vehicle  and  to  move  with  it,  the  reference  frame 
is  rotating  and  the  equations  become 
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F  =  m  dtf/dt  +  m  (D  X  \?c  FTT-3 

M  =  ah/at  +  u  X  h  TTT-A 

Resolution  of  each  of  these  total  vectors  into  its  components, a  long  and 
about  the  axes,  gives  the  Eulerian  equations  of  motion  as: 


Fx  =  m(C  +  QW  -  RV) 
Fy  =  m(V  +  RU  -  PW) 
Fz  -  m(W  +  PV  -  gU) 


L 

=  P'x 

-  Rlxz  + 

QR(i2  - 

ly) 

1 

3 

X 

N 

M 

^ly 

+  PR  ( I  x  • 

-iz)  -  R2 

1  xz 

+  p2 1 
+  r  ‘xz 

N 

=  Rlz 

-  Plxz  + 

3 

*< 

1 

lx) 

+  QRixz 

where  Fx,  F  ,  and  Fz  are  the  externally  applied  forces  along  the  x,  y,  z 
axes;  L,  M,  and  N  are  the  externally  applied  moments  about  the  x,  y,  and 
2  axes;  U,  V,  and  W  are  vehicle  velocities  along  the  x,  y,  and  z  axes;  P, 

Q,  and  R  are  rotational  velocities  about  the  x,  y,  and  2  axes;  !x,  ly, 
lz  are  moments  of  inertia  about  the  respective  axes;  lxz  is  the  product 
of  inertia  about  the  y-axis,  and  the  xz  plane  is  taken  as  a  plane  of  sym¬ 
metry.  Equat ion s ILL- 5  are  then  the  general  non-linear  equations  of 
mot  ion. 

Li  near i zat i on 

The  equations  of  motion  may  be  linearized  by  the  assumption  that 
all  perturbations  away  from  the  equilibrium  state  will  be  small..  The  velo¬ 
cities  of  Equations  XEI-5  may  then  be  taken  as  summations  of  the  steady 
state  values,  indicated  by  the  subscript  or  plus  the  perturbation  quanti¬ 
ties,  indicated  by  lower  case  letters. 
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U  -  Uo  +  u 
V  =  Vo  +  v 
W  =  wQ  +  w 

m-6 

p  =  pQ  +  p 
Q  =  Qo  +  q 


R  -  Rq  +■  r 

Substitution  of  Equations  LJ_L_6  into  Equations  I  TIL-5,  and  the  neglection 
of  products  and  squares  of  perturbation  quantities  gives,  with  the  assump¬ 
tion  of  wing  level  steady  flight  with  all  velocity  components  zero  except 
Uo  and  W0  ,  : 

Fx  =  m[u  +  W0q] 

Fy  =  m[v  +  U0r  -  WQp] 

Fz  =  m[w  -  l^q] 

Ed- 7 

L  =  P  i  x  "  ^  I  xz 

M  =  qly 

N  *  rlz  -  plxz 

The  Euler  axes  are  assumed  to  rotate  from  one  position  to  another  in  the 
order  0,  cp  as  shown  in  Figure  II.  The  relation  between  the  angular 
velocities  and  these  Euler  angles,  which  are  small  since  the  perturbations 
are  sma  I  I ,  are  the  r- 


p  =  cp 

q  =  9  m-8 


r  =  * 

Gravity  and  Aerodynamic  Forces 

The  gravity  force  on  a  body  in  the  atmosphere  acts  along  a  line 
through  the  center  of  the  earth.  If  the  equilibrium  Euler  axes  are 


7 
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disp laced  form  this  line  by  the  angles  0O  and<pt>,  then  the  steady  flight 
components  of  the  weight  force  are  given  by 

9xo  =  W  sin  0O 

9y0  =  °  xn-9 

gz^  =  W  cos  0O  cos  Cfo 

where  the  wing  level  steady  flight  condition  <po  =  0  has  been  i ncorporated. 
In  disturbed  flight,  these  components  must  be  rotated  through  the  Euler 
angles  t,  0,  cp.  This  gives  the  weight  components  along  the  disturbed  axes, 
with  the  disturbance  assumed  small,  as: 

g*  =  -  V/  sin  8C  -  W  cos  0O  0 

gy  —  W  sin  0O  i|f  +  W  cos  80  cp  I  IT- 1 0 

g2  =  -  W  sin  G0  0  +  W  cos  0O 

The  aerodynamic  forces  and  moments  acting  on  the  vehicle  may  be 

expanded  in  a  Taylor's  series  about  an  equilibrium  position  as  follows 
G  =  G0  +  (Ga)g  a  +  (Gjjj)0  8  +  (Gu)0  u  +  ...  +  H.O.T. 
where  (Ga)0  =  dG/da  evaluated  at  the  equilibrium  conditions,  etc.  and  H. 
O.T.  indicates  higher  order  terms.  The  a,  0,  u,  etc.  are  the  parameters 
with  which  the  aerodynamic  forces  and  moments  are  known  to  vary.  They  are 
perturbation  quantities  and  the  higher  order  terms  which  consist  of  powers 
of  the  quantities  of  order  two  or  greater  may  be  neglected.  If  X,  Y,  and 
Z  are  used  to  indicate  the  aerodynamic  forces  along  the  x,  y,  and  z  axes 
and  L,  M,  and  N  are  used  to  indicate  the  moments  about  the  axes,  the  expan¬ 
sions  of  the  aerodynamic  forces  and  moments  in  terms  of  parameters  with 

which  each  is  known  to  vary  in  quasi -steady  flow  are  then 

X  =  XQ  +  Xuu  +  Xqq  +  Xww  +  X^w  +  X*u  +  X^q  HI- 1  la 

y  -  Y°  +  v  +  V  +  V  +  V  +  Yv*  +  V 


m-i  ib 
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Z  =  ZQ  +  Zuu  +  Zqq  +  Zww  +  Z^u  +  +  Z^w  ITU—  1  lc 
L  -  Lq  +  Lrr  +  Lvv  +  Lpp  +  L^r  +  L-v  +  L*p  JXL-I Id 
M  =  1^,  +  f^u  +  Mqq  +  Mww  +  M^Ci  +  M^q  +  M^w  HI- 1  le 
N  =  Nq  +  Nyv  +  N«v  +  Nrr  +  N*f  +  Npp  +  N*p  +  Nww  HI- 1  If 


Then,  if  the  Euler: on  axis  system  is  specified  as  a  stability  axis  system, 
WD  =  0  and  all  partial  derivatives  with  respect  to  rates  of  change  of  velo¬ 
cities  are  zero.  The  equations  of  motion  obtained  from  Equations  lll-ll, 

10  and  7  become: 

mu  =  X0  +  Xuu  +  Xqq  +  Xww  -  W  sin  y0  -  W  cos  yo0  +  f xz  +  TXq 
m(v  +  U0r)  =  Y0  +  Yrr  +  Yvv  +  Ypp  +  W  sin  y0iJi'+W  cos  yDcp  +  Tyz 

m(w  -  L^q)  =  Z0  +  Z^u  +  Zqq  +  Zww  -  W  sin  y08  +  W  cos  yQ  +  Tzz  +  T2q 

P!x  “  r  I  xz  =  Lo  +  Lrr  +  Lyv  +  Lpp  +  TL ^  EH- 1 2 

qly  =  Nb  +  Muu  +  Mqq  +  Mww  +  f„z  +  TMq 
r I z  -  plxz  =  Nq  +  Nyv  +  Nrr  +  Npp  +  Nww  +  fNz 
where  the  Tq  terms  represent  variable  force  parameters  and  the  Tq  steady 
state  cable  force  values. 

For  the  system  in  equilibrium  flight,  a  I  I  of  the  perturbation 
terms  and  time  derivatives  are  zero.  The  equi librium  equations  are  then 

Xo  -  W  sin  Yo  +  =  0 

Z0  +  W  cos  y0  +  Tz  =0 

0  m-13 

**>  +  tHq  -  0 

Y0  =  Lo  =  No  =  0 

Substitution  of  Equations  IU-13  into  Equations  I  IX-12;  division  of  the 
first  three  equations  of  Equation  HI-12  by  m,  and  division  of  the  roll 
equation  by  Ix,  the  pitch  equation  by  ly  and  the  yaw  equation  by  lz;  and 
substitution  of  Equations  I  I  I -fl  yields  the  final  form  of  the  equations  of 


mot i on . 
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u  -  Xuu  -  Xq0  -  Xww  +  g9  cos  yQ  -  Txz  =  0 
w  -  Uq 9  -  ZjjU  -  Zq5  -  Zww  +  g9  sin  Yo  -  Tzz  =  0  HI- 1 4 
8  -  Muu  -  Mq§  -  M„w  -  T„z  =0 

v  +  'Jot  -  Yn|  -  Yyv  -  Ypcp  -  t9  sin  y  -  gcp  cos  y0  -  Ty  =  0 
‘tp  -  V  Uxz/Ix)  -  LrJ  -  Lvv  -  Lpcf>  -  Tl  =0  m-15 

t  -  cp  (lxz/l2)  -  Nrt  -  Nvv  -  Npcp  -  Nww  -  =  0 

The  parameters  in  Equations  IEI-I4  and  I  I  I  - 1 5  obtained  by  dividing  the 
partial  derivatives  of  the  aerodynamic  forces  by  inertial  terms  (i.e.  Xy  = 
I/m  X„;  =  I / ! y  M(J;  Lv  =  l/lx  Lr;  Nr  =  l/lz  Nr;  etc.)  are  termed  dimen¬ 

sional  stability  derivatives. 

The  aerodynamic  forces  and  moments  may  be  expressed  in  the  form 

F  =  Cp  1/2  pV3  S 

and  M  =  C{>.  1/2  PV2  S  C 

where  Cp  and  CM  are  dimensionless  coefficients,  p  is  the  fluid  density, 

V  is  the  total  vehicle  velocity  and  S  and  C  are  characteristic  area  and 
length.  Utilizing  these  definitions  in  obtaining  the  partial  derivatives, 
the  dimensional  stability  derivatives  of  Equations  111-14  and  111-15  may 
be  obtained.  A  complete  derivation  of  these  terms  may  be  found  in  Re¬ 
ference  10.  The  dimensional  stability  derivatives  of  interest  here  are 
presented  in  Table  I. 

Cable  Forces 

The  transverse  vibrations  of  the  cable  were  seen  in  Section  H, 
for  the  case  of  time  varying  boundaries,  to  vary  with  both  time  and  posi¬ 
tion  along  the  cable.  The  cable  form  can  then  be  considered  to  consist 
of  a  traveling  wave  superposed  upon  a  standing  wave.  In  particular,  since 
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TABLE  I 

DIMENSIONAL  STABILITY  DERIVATIVES 


Lonqi tudi na i 

La+era 1 

xy  =  (pus/fnx-  c0u  -  C0) 

Yv 

=  (pUS  /2m)  Cyg 

Xw  =  (PUS/2m) {CL  -  C0  ) 

Yr 

=  (pUSb/4m)  Cyr 

Xq  fix  0 

YP 

=  (pUSb/4m)  Cyp 

Zu  =  CpUS/m) (Ct  -  CL) 

Lv 

=  CpUSb/2 Ix)  C£& 

Zw  =  {pUS/2m)(-  CLa  -  C0) 

^■r 

=  CpU$b*/4lx)  C£ 

2q  =  CpUSC/4m)  CLq 

Lp 

=  CpUSbs/4Ix) 

Mu  =  (pUSc/Iy)(Cmn  +  Cm) 

Nv 

*  (pUSb/2tz)  Cng 

Mw  =  (pUSc/2ly)  C% 

Nr 

=  CpUSbs/4lz)  C0( 

Mq  =  (pUScs/4ly)  Cpr^ 

Nw 

=  (pUSb/21z)  Cn„ 

•A 
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over  80  percent  of  the  energy  in  a  fixed  end  vibrating  cable  is  contained 
in  the  fundamental  mode,  the  cable  form  may  be  considered  as  a  traveling 
wave  superposed  on  the  fundamental  which  in  effect  is  the  equilibrium  cable 
shape.  In  addition,  the  traveling  wave  has  been  shown  to  be  damped  in 
both  the  upstream  and  downstream  direction  at  low  velocities  which  are  con¬ 
sidered  in  this  report.  Consequently,  the  traveling  wave  may  be  neglected 
and  the  cable  may  be  assumed  to  vary  from  one  equilibrium  shape  to  another 
with  perturbations  in  the  towed  vehicle  position. 

Therefore,  to  determine  the  cable  contributions  of  the  dynamic 
stability  of  the  system,  the  following  assumptions  are  made: 

1)  The  cable  shape  is  that  of  its  fundamental  mode 
and  remains  that  way  through  out  all  reasonable 
perturbations  from  the  equilibrium. 

2)  Small  angle  approximations  are  valid. 

The  cable  tension  T  may  be  resolved  into  the  horizontal  and  ver¬ 
tical  components: 


T  cos  X 


HI-16 


and 


T  sin  X  HI-17 

respectively.  The  equations  of  motion  are  such  that  the  variation  of 

these  tension  components  for  a  given  disturbance  must  be  knownfi.e. 

dT  cos  X/dz,  etc.  must  be  determined).  Differentiation  of  Equation  HI- 

16  with  respect  to  z  gives: 

d  T  cos  X  =  -  T0  si  n  X0  di.  +  cos  X  dT 
dz  dz  dz 

dT/dz  and  dT/dX  are  determined  as  follows.  The  derivative,  dT/dX,  is  ob¬ 
tained  from  Figure  8  by  measuring  the  slope  of  the  point  corresponding  to 
the  known  values  of  tp,  the  equilibrium  cable  angle,  a,  the  cable  parameter 
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and  the  ratio  T0/T.  The  expression  dX/dTj  may  be  obtained  from  Equation 
IX- 10 1  by  first  replacing  cp  by  the  i  dent  i  f  i  cat  i  on  cp  =  90  -  X  and  after 
simple  algebraic  manipulation,  differentiation  with  respect  to  T|.  Multi¬ 
plication  of  dT/dX  by  dX/dT)  yields  dT / dT)  which  is  the  non-dimensional  form 
of  dT/dz. 

It  is  noted,  hcwever,  that  this  development  did  not  account  for 
any  changes  in  cable  contributions  with  a  change  in  length  of  the  cable. 
This  can  Ibe  seen  from  Figure  7;  for  although  the  cable  becomes  shorter, 
the  attachment  point  remains  fixed  and  thus  the  value  of  the  cable  deri¬ 
vatives  remain  constant.  Therefore,  for  a  given  displacement  Az,  the 
change  in  X  is  approximately  a  function  of  i  as  follows  (see  Figure  12): 

dX  sec  X  .  cos  X  sa  X  HI-  18 

dz  l  t 

This  effect  was  computed  for  a  range  of  a's  of  2°  through  20°  and  a  number 
of  rope  lengths.  The  results  of  this  addition  to  the  cfcp/dz  term  previously 
developed  are  discussed  in  a  later  section  and  are  most  significant. 

The  cable  contributions  in  the  longitudinal  equations.  Equations 
ill- 1 4,  with  the  perturbation  angle  a  neglected  in  comparison  with  the 
large  angle  X  are: 

Tx  =  l/m[cos  X0  (3T/3z)  -  T0  sin  X0OX/dz)] 

Tz  =  l/m[sin  X0 (dT/dz)  +  T0  cos  XoCdX/dz)]  HI-19 

Tm  =  l/ly[e(T0  cos  Xc{dX/dz)  +  sin  X0OT/dz))  -  a(cos  X0  (dT/dz) 

-To  sin  x0 cax/az)) 

The  geometric  parameters  b  and  a  may  be  determined  from  Figure  13. 

The  terms  TY ,  TL ,  and  T^  of  Equations  111-15  could  be  determined 
from  a  consi deration  of  the  variation  of  the  cable  shape  that  results  from 
a  lateral  displacement  in  the  same  manner  as  the  longitudinal  case  was 
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developed.  However,  the  contributions  due  to  this  variation  in  cable  shape 
is  small,  since  the  displacements  to  be  considered  are  small,  and  in  ad¬ 
dition  they  are  stabilizing  contributions.  Consequently,  the  effect  of 
neglecting  changes  in  cable  shape  will  be  minor  and  conservative. 

The  cable  is  therefore  assumed  to  remain  straight  in  the  lateral 
directional  planes.  The  lateral  directional  cable  forces  may  then  be  de¬ 
termined  from  the  geometry  of  Figures  14  and  15,  with  all  angles  assumed 
to  be  sma 11,  to  be 

TY  =  l/m[(l  +  e/j t)3  -  y*/i  -  (I  +  a/n)cp][T0  +  (ST/Bz)z] 

Tl  =  a/lx[(l  +  e/i)p  -  y */ 1  -  (I  +  a/n)cp][Tc  +  OT/9z)z]  m-20 

Tn  =  b/lx[CI  +  e/£)g  -  y */£  -  Cl  +  a/n)cp][T0  +  OT/dz)z] 

The  complete  equations  of  motion  may  then  be  obtained  when  necessary  by 

incorporating  Equations  I  1 |- 1 9  and  20  into  Equations  1 11-14  and  15  respec¬ 
tively.  There  are  then  six  equations  with  eight  time  dependent  variables. 
This  is,  of  course,  an  undesirable  situation.  |t  may  be  alleviated,  how¬ 
ever,  by  considering  the  rate  of  side  displacement  and  the  rate  of  climb 
relations.  The  rate  of  climb  z  is  given  by 

z  =  U  sin  y 

where  y  is  the  flight  path  angle.  The  small  perturbations  assumption  and 
the  assumption  that  y0  is  very  small  or  zero  leads  to 

z  =  Uq  y 

Integration  and  the  substitution  of  y  =  9  -  a  then  gives  the  desired  form 

z  =  j  UqCS  -  a)  dr  IH-21 

fo 

Likewise,  the  rate  of  side  displacement  is  given  by 


y*  =  u0  * 
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where  y*  is  the  displacement  perpi ndi cu lar  to  the  equi librium  x-axi 
Upon  integration,  this  gives 

rt 

y*  =  P  U0  t  dT 

t0 

With  the  addition  of  Equations  1 1 1 -2 1  and  22,  the  equations  are  in 
plete  form  which  may  be  solved. 


s. 


HI -22 


a  com- 
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I?  LONGITUDINAL  DYNAMICS 

The  study  of  the  longitudinal  dynamics  of  a  towed  vehicle  pro¬ 
ceeded  in  the  normal  manner  from  simplified  analytical  studies  to  the  more 
complicated  studies  which  required  the  use  of  an  analog  computer.  The 
simplified  studies  were  designed  to  indicate  important  vehicle  parameters 
and  to  determine  the  value  of  different  degrees  of  approximations  of  the 
cable  tension  variation. 

In  order  to  establish  relative  magnitude  of  the  various  parameters, 
it  was  necessary  to  specialize  the  studies  to  a  particular  basic  configu¬ 
ration.  The  basic  configuration  chosen  is  shown  in  Figure  13.  The  static 
stability  derivatives  for  this  configuration  were  obtained  from  preliminary 
wind  tunnel  dafa  of  tests  performed  at  the  David  Taylor  Model  Basin.  The 
equilibrium  flight  parameters  were  calculated  from  Equations  HT-13.  The 
procedure  followed  is  described  in  Appendix  B  and  the  parameters  for  sev¬ 
eral  conditions  are  listed.  The  basic  flight  condition  of  170  m.p.h.  at 
sea  level  and  a  cable  length  of  250  feet  were  chosen  for  the  majority  of 
the  studies. 

The  results  are  presented  in  the  form  of  stabi  lity  boundaries 
where  applicable,  and  the  effects  of  variations  in  the  other  parameters 
are  di scussed. 

Two  Degree  of  Freedom  Studies 

The  two  degree  of  freedom  studies  were  performed  with  the  assump¬ 
tions  that  both  the  tension  in  the  cable  and  the  angle  between  the  cable 
and  the  body  centerline  (X)  remain  constant.  The  cable  contributions  of 
Equations  l l 1  - 1 4  are  then 
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Txz  =  (-  T0/m  sin  Xo)a  =  -  Tx* 

Tzz  -  CT0/m  cos  X0)a  =  Tz* 

Tp-Z  =  T0/ly(e  cos  X0  +  a  sin  X0)a  =  Tm* 

With  the  flight  path  angle  y0  taken  as  zero,  and  w  =  U0a,  the 
Laplace  transformed  equations  of  longitudinal  motion  become: 

XqS0  +  (Xu  —  s) a  +  (Xa  -  Tx*)a  —  0 
<Zq  +  U0)se  +  Zuu  -  CUbs  -  Za  +  Tz*)a  =  0  U-l 

-(s2  -  MqS) 9  +  Muu  +  Ma  +  T„*  =  0 
The  two  fundamental  modes  of  motion  of  most  vehicles  of  this 
type  are  the  short  period  and  hne  phugoid.  The  short  period  is  highly 
damped  and  the  velocity  of  the  vehicle  changes  very  little  before  its  ef¬ 
fect  becomes  negligible.  Therefore,  the  short  period  mode  may  be  studied 
approximately  by  setting  u  =  0  and  studying  the  lift  and  moment  equations 
Setting  y0  =  0,  the  character i sf i c  equation  for  the  short  period  mode  is 
then : 

(Za+U0)s  -  UoS2+  (Z  -  Tz*) 

^  =  o  nr-2 

-(s2  -Mqs)  ^  +  Tm* 

Expanding,  this  becomes: 

Css  -[Mq  +  (Za/Uk,>]~  tJ/uJs  +  (Za  -  Tz»)  -  (Zq/U„)  (Hj+r,*) 

-  -  T„*]  =  0  EE-3 

Equation  ET-3  is  of  the  form  S2+  2 C,  u)rS  +  cun2  =  0,  and  we  know  that  for 
stabi  liiy  2£  c-Jn  =  -[Mq  +  Za/Uo  -  T^/Uq]  must  be  positive.  Therefore,  the 
criterion  for  stabi I i ty  of  this  approximation  to  the  short  period  mode  is 

[Mq  +  (Za/U0)  -  <Tz*/Uo>  <  0  ET-4 

An  order  of  magnitude  investigation  of  the  parameters  involved  indicates 
that  Mq  <  0,  Zg/Ur,  <  0  and  Tz*/Uo  >  0.  Therefore,  [Mq  +  Za/U0  -  Tz*/Uol 
<  0  and  the  approximation  to  the  short  period  oscillation  is  stable. 
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The  motion  in  the  Phugoid  mode  consists  mainly  of  altitude  and 
speed  changes  with  very  small  a  variations.  This  mode  may  then  be  studied 
by  setting  a  =  0.  In  addition,  the  motion  is  slow  so  that  the  inertia 
terms  are  negligible  and  since  Mq  and  M,j  are  usually  very  small,  the  pitch¬ 
ing  moment  equation  may  be  neglected.  The  characteristic  equation  of  the 
Phugoid  mode  is  then: 

X  s  -  g  -  s 

=  0  I?-5 

(Zq  +  Uc)s  Zu 

where  again  y0  =  0.  Expansion  then  yields: 

s3  +  (  ZUXQ  -  Xu)s  -  qZu  =  0  H-6 

Uo+Zq  Uo  +  Zq 

The  criterion  for  Phugoid  stability  is  then  that 

[(Z^q/Uo+Zq)  -  Xj  >  0  IY-7 

An  order  of  magnitude  investigation  indicates: 

Xu  <  0  Xq  <  0 

Zu  <  0  Zq  <  0 

U0  +  Zq  >  0 

Therefore,  [ZyXq/CUo  +  Zq)  -  Xu]  >  0  and  the  approximation  to  the  phugoid 
mot  ion  is  stab  I e . 

Three  Degree  of  Freedom  Analytical  Studies 

The  dependence  of  the  phugoid  approximation  upon  the  parameters 
Mq  =  0  and  Mu  =  0  makes  it  desirable  to  check  the  complete  three-degree 
of  freedom  stabi lity  for  this  configuration.  To  check  the  importance  of 
the  q  derivatives,  a  typical  flight  condition  of  170  mph  and  an  aerodynamic 
center  location  of  2.5  feet  aft  of  the  center  of  gravity  have  been  used. 

A  one  degree  change  in  the  effective  angle  of  attack  due  to  a  positive  q 
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rotation  about  the  c.g.  requires  the  creation  of  a  vertical  velocity  of 
4.37  feet/sec.  at  the  aerodynamic  center.  The  vertical  velocity  at  the 
aerodynamic  center  due  to  q  is  given  by 

w  -  q ^xac  -  xcq^ 

Solving  this  equation  for  q,  it  is  seen  that  a  pitch  velocity  of  100  deg/ 
sec  is  required  to  give  a  one  degree  change  in  effective  angle  of  attack. 
Therefore,  if  may  ho  assumed  that  CL^  =  CD^  =  CH(^  =  0,  and  correspondingly, 

Zq  '  ><q  -  Mq  =  °- 

The  dimensional  stability  derivative  My  is  given  by: 

Mu  =  (pU0Sa/lyy)(CMu  +CMo) 

The  coefficient  CM  is  a  Mach  number  effect  and  may  of  course  be  assumed 
zero  in  the  speed  range  of  interest.  CH  is  the  aerodynamic  coefficient 
of  pitching  moment  in  the  equilibrium  flight  condition  and  may  be  suspec¬ 
ted  to  have  an  appreciable  negative  value  to  balance  the  moment  created 
by  the  line  tension.  The  value  of  CMq  is  obtained  from  the  appendix  as 
CHq  =  -.28  which  as  anticipated  is  an  appreciable  negative  value  and  the 
three  degree  of  freedom  stabi lity  must  be  checked. 


Constant  Tension  -  Constant  X 

The  three  degree  of  freedom  characteristic  equations  correspond¬ 
ing  to  the  assumption  of  constant  tension  and  constant  X  is  obtained  from 
Equations  TV- 1  as: 

Xqs  “  s  -  'x* 

CZq  +  (J0 )  s  Z„  “(Uo  s  -  Za  +  Tz*)  =  0  H-8 

-(ss4WqS)  Mu  Ha+T„* 

Expanding  the  determinant  yields  the  equation  in  the  form 

As4  +  Bs3  +  Cs2  +  Ds  +  E  =0  H-9 
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where: 

A  =  I 

B  =  -[X.j  +-  (Za/H>)  -  {T2*/Uc  )  ] 

C  =  [XuCZa/Uo  -  T2*/Uo)  “  ^  -  Zu(Xa/Uo  -  Tx*/Uo)  -  TH*  ] 

B  =  [MuCg  -  Xa  +  Tx*)  +  Xu(Ma  +  T H*) ] 

E  =  gtZ^Ma/Uo  +  TH*/Uo)  -  M!u(Za/Uo  -  Tz*/Ub ) ] 

With  the  equation  in  this  form,  the  stabi  lity  may  be  checked  by  applying 
Routh's  stability  criteria  for  a  quartic.  Routh’s  criteria  tor  a  quartic 
state  that  for  stabi  lity,  the  constants  of  the  equation  must  meet  the  con 
ditions  that 

A  >  0,  B  >  0,  D  >  0,  E  >  0 

and 

D(BC  -  AD)  -  B2E  >  0 

To  apply  these  criteria  to  equation  DT-9,  the  dimensional  stability  deri¬ 
vatives  must  be  evaluated  numerically.  For  the  equilibrium  conditions 
stated  previously,  these  values  are: 

X„  =  -  0. 197;  Zu  -  -  0.094 

Mu  =  -  0.052;  Xa/Uo  =  +  0.047 

Za/Uo  =  -  L9V8;  N^/U0  =  -  0.78 

Tx/Uo  =  +  0-0818:  Tz/U0  =  +  0.098 

T„/U0  =  +  0-0456 

The  coefficients  of  Equation  H-9  may  be  evaluated  as: 

A  =  +  !  >0 

B  =  +  2.28  >  0 
C  =  +  129  >0 
D  =  +  22  >0 

E  =  -  1.39  <  0 


-52- 


«■ 


and  it  may  be  seen  that  the  first  part  of  Routhfs  criteria  is  not  satisfied 
for  £  <  0.  Checking  the  second  part, 

0  (BC  -  AD)  -  B3E  =  +  7577  >  0 
and  the  criterion  is  satisfied. 

One  of  Routhfs  criteria  is  not  satisfied  and  the  system  is  there¬ 
fore  unstable.  However,  the  second  criterion  is  satisfied  and  it  may  there¬ 
fore  be  suspected  that  the  instab i lity  is  weak.  This  type  of  instability 
is  characterized  by  the  depletion  of  the  phugoid  mode  from  an  oscillation 
into  two. real  roots,  one  of  which  is  convergent  and  the  other  divergent. 

This  phenomenon  is  here  directly  due  to  the  strong  Mu  term.  Ordinari ly, 
in  piloted  or  otherwi ze .contro I  led  aircraft,  this  divergence  resulting 
from  the  Mu  term  would  cause  no  great  concern,  for  it  is  normally  a  slow 
divergence  and  can  readi ly  be  controlled.  In  this  case,  however,  control 
is  not  possible  and  the  effects  of  further  variations  in  the  cable  contri¬ 
bution  must  be  considered.  If  necessary,  this  effect  may  be  alleviated  by 
shifting  the  center  of  gravity  forward  or  by  otherwize  changing  the  con¬ 
figuration  so  that  the  equilibrium  position  moment  due  to  line  tension 
i s  reduced. 

Tension  Constant,  X  Semi -constant 

The  next  step  toward  the  development  of  the  actual  system  anal¬ 
ysis  is  to  again  consider  that  T  remains  constant  at  its  equi librium  value, 
but  that  the  angle  X  is  semi -constant.  What  is  meant  by  considering  X  as 
a  semi -con  <=Tant  can  best  be  seen  by  considering  the  vehicle  at  its  initial 
displaced  position  following  a  disturbance.  According  to  the  basic  no 
line  harmonics  assumption  for  these  sfudies  mentioned  previously,  the  line 
would  row  assume  a  tension  and  an  angle  X  with  the  relative  wind  that  would 
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agree  with  some  equilibrium  vehicle  configuration  at  that  position  in  space. 
The  tension  is  now  assumed  constant  at  the  original  equilibrium  value  T0 . 
The  angle  X  at  this  initial  displaced  position  now  has  some  value  \0  +  . 

Then  since  X  =  a  +  K,  Xi  =  ax  +  X x .  The  variation  in  X  from  Xx  is  then 
assumed  to  consist  only  of  a  variations,  and  X  remains  constant  at  X  =  X* , 
such  that  X  at  any  later  time  is  given  by  X  =  a  +  or  X  =  a  +  X  where 
X  is  considered  constant.  This  then  allows  transfer  functions  +o  be  estab¬ 
lished  in  much  the  same  way  as  control  deflection  transfer  functions  are 
estab 1 i shed. 

The  equations  of  motion  with  this  assumption  and  =  Zq  =  Mq  = 

Yo  =  0  are  then; 

-  g0  +  XyU  -  u  +  Xaa  -  T^a  =  Tx*  X 

U09  +  Zuu  +  Zaa  -  lU  -  Tz*  a  =  Tz*  X  EMO 

6  -  Muu  -  l^a  -  T„*  a  =  TH*  X 

Laplace  transforrni  ng,  these  equations  become: 

-g  Xu  -  s  xa  Tx*  ®  Tx 

Uos  Zu  Za  -  Tz*  -  UoS  u  Tz  X  DT-I  I 

sa  -My  -<Ma  +  TM*)  a  T„ 

The  transfer  functions  of  interest  may  now  be  defined  as: 

9j^I  =  No;  Llisi  =  %;  ssJUi  = 

X(s)  A  X(s)  A  XCs)  A 

The  denominator  is  then; 

A  =  U0  s4  -[(Za  -  Tz*/Uo)  +  Xy]  s3  +[Xu(Za  -  Tz*/Ub)  -  Zu (Xa-Tx*/Uo ) 

“  Ma  -  TM*]  ss  +  [XUC^  +  T„*)  +  Mu(g  -  Xa  +  Tx*)  ]  s  +  g[Zu(^  +  T**/^) 

+  Mu(Za  -  Tz*/U0)]  W-12 
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No  =  U0  T„»  s2  -[  (TX*/TH)MU  +  (Tz*/TH)(Ma/U0)  +  Xu  +  (Za/4>)]  s 
+  Tx*/TH(MuCZa/Uo)  -Zu(^/Uo>)  +  TZ*/TM  (Xu  (Vuo>  "  MuCXa/Uo)) 

+  X.j(Za/iJ0)  -X.,(Xa/Jo) 

Nu  =  -  L'0  Tx*  s3  +  [az*/Tx*ma/Ub)  -  (Za/Uo)]  s2  +  [CTH*/TX*)  • 

(g  -  Xa)  -  ^]s  -  g[(T,/Tx)(l^/Uo)  +  (TH/TX)  CZa/Uo)]  H-13 

Na  =  “  Tz  ~  &u  -  CTx/Tz)Zu  +  Lb  Tj  s2  +  [Ub  CTX/TZ)MU 

+  UoaH/Tz)Xu]s  +  g[M.j  +  ZuCTm/T2)] 

Substituting  the  values  of  the  dimensional  stability  derivatives  presented 

previously,  the  transfer  functions  are: 

6(5)  =  _ (11.4)  (s  ±  .209) (s  +  3.735) 

XCs)  Cs  +  .2185) (s  -  .0315) Cs2+  2.083  s  +  183.65) 

iii^L  =  -  20.45(s  +  .318) (s2  +  1.717  s  +  206.454) 

X  Cs)  (s  +  .2185)  (s  -  .0315)  (s2  +  2.083  s  +  183.65)  -LV-‘14 

aSs l  =  -  0.0982(s  -  2849.89) Cs2  +  .0119  s  +  .0011) 

XCs)  (s  +  .2183) (s  -  .0313) (s2  +2.003  s  +  183.65) 

Approximate  zero-pole,  s-p lane  locations  of  these  transfer  functions  are 

shown  ir,  Figure  16.  The  system  will,  of  course,  operate  at  points  determined 

by  the  gain  of  the  transfer  functions.  The  gain  in  each  case  is  seen, 

from  Equations  DT-13,  to  be  a  direct  function  of  the  line  tension.  In 

view  or  this  dependence  or.  line  tension  and  the  fact  that  allowing  X  to 

be  semi -constant  does  not  alter  the  location  of  the  pole  in  the  right  half 

plane  appreciably,  direct  attention  W3S  focused  on  the  variable  tension 

case. 


Static  Effect  of  Variable  Tension 

Essentially,  the  tension  in  the  line  will  be  expressed  by  an 
expansion  about  the  equi librium  tension.  This  expression  has  the  basic 
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t  =  t0  +  di/ax  x  nr- i 5 

Some  insight  as  to  the  expected  results  may,  however,  be  obtained  by  look¬ 
ing  at  the  static  stability  contribution  of  the  dT/dX  term  which  is  known 
to  be  negative. 

Consider  a  situation  where  the  vehicle  is  disturbed  by  an  ini¬ 
tial  change  in  angle  of  attack  a  at  some  time  tx  .  Let  the  subscript  0 i n- 
dicate  equilibrium  conditions  and  the  subscript  :  the  conditions  at  time 
tT .  Referring  to  Figure  17,  the  moment  about  the  c.g.  due  to  the  line 
tension  is  positive  and  is  given  by: 

M,  =  a  T  sin  g  I?-  1 6 

Since  fir  is  positive,  the  condition  for  static  stability  is: 

MTl  <  MTo  H-17 

The  change  in  position  of  the  vehicle  that  results  from  the  change  in  a 
is  not  being  considered  (this  would  be  a  dynamics  study),  only  the  condi- 


tions  at  the  time  the  change 

in  a  is  applied.  Therefore, 

X^  —  \q  and 

since  X  =  g  +  £,  where  §  is  £ 

1  constant  dependent  only  on 

the  geometry  of 

the  vehicle. 

sin  gi  =  si n  g0 

ft 

1 

00 

From  equation  TV- 16.  we  have 

for  and  MTq 

Mtj  =  Tj  a  si  n  gy 

IT- 19 

MTq  -  T<j  a  sin  g0 

U-20 

Then  di vi di ng 

=  T,  a  sin  g,  =  Ty 

DT-21 

Mr 

T  0 

T0  asingg  T0 

Now  from  the  basic  relation. 

equation  U-15: 

T 
*  1 

=  T0  h  9T/SX  X 
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Now  X  =  Aj  -  A0 

and  A1  =  \1  +  ax  =  Xx  +  do  +  a 

Aj  =  Xo  +  ao 

Then  substituting  Equations  TV-? 5  and  24  in  Equation  IT-22 

X  —  Xj  ^  cl  —  Xq  , 

ar.d  si  nee  Xx  =  X0 

X  =  a  >  0 

Then  since  9T/3X  <  0 


aT/dX  X  <  0 

and  T0  +  aT/dX  <  T0 

Therefore,  from  Equations  TV-  I  5  and  IS- 27,  it  is  seen  that 

Tj.  <  To 

and  from  Equations  H-2 I  and  IT-28, 

M^/M^  =  Ti/T^  1 

Consequently,  relation  DZ-17 


IV-22 

IY-23 

nr-24 


IV-25 

IT-26 

IE-27 

ESr-28 


«rl  <  Mt0 

is  satisfied  and  the  line  tension  contribution  is  statically  stabilizing. 
It  appears,  therefore,  at  least  from  static  stability  considerations,  the 
inclusion  of  variable  line  tension  will  stabilize  the  system. 


Complete  Three  Degree  of  Freedom  Studies 

The  foregoing  analyses  have  indicated  that  the  basic  system  may 
be  expected  to  be  dynamically  stable  when  the  complete  effect  of  the  cable 
forces  are  included  in  the  analysis.  This  may  of  course  be  done  analytic¬ 
ally;  however,  the  ultimate  aim  is  a  study  of  the  variations  of  several 
of  the  parameters  away  from  their  basic  values.  Therefore,  an  analog  com¬ 
puter  was  utilized  for  the  remainder  of  the  study. 
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The  longitudinal  equations  of  Equations  I  I  1-14  may  be  put  in  a 
form  better  suited  to  a  parametric  study  on  an  analog  computer.  This  may 
be  accomplished  by  substitutions  of  the  non-dimensional  stability  deriva¬ 
tive  forms  of  Table  I  for  the  dimensional  stability  derivatives  in  the 
equations  and  defining  the  operator  D  =  d/dt.  The  Equations  El  1  - 1 4  then 
become 

[D  +  CpUS/m)C0„]u  -  C0U2S/2m(a  -  C^la  +  g8  -  Txz  =  0 
(pUS/m)CLu  +  [U0D  +  (pUS/2m)(CLa  +  CDo)]a  -  Dl^e  -  Tzz  =0  IT-29 

(pUSc/ly)CHo  -  (PU3SC/2ly)CMa  a  -  [D2  -  CDUSc3/4| y)CMq  D]0 

-  THz  =  0 

Equations  IT-29  along  with  Equations  HI-19  and  1 1 1 -2 1  were  solved  on  the 
analog  computer  to  obtain  the  results  discussed  i  r.  the  remainder  of  this 
section. 

All  of  the  aerodynamic  parameters  of  Equations  IT-29  were  varied 

to  some  degree.  The  three  parameters  of  primary  importance  were  determined 

to  be  C„a,  CHq  and  Cw„.  The  first  two,  C„a  and  C„q  are  normally  important 

in  common  airplane  dynamics  and  the  effect  was  as  suspected.  The  parameter 

C„o  is  somewhat  peculiar  to  towed  vehicles  since  it  has  an  unusually  high 

negative  value  as  discussed  in  an  earlier  section. 

The  changes  in  the  other  parameters  in  general  have  very  little 

individual  effect.  In  particular,  any  variation  that  results  from  either 

C5o  or  CUQ  changes  can  be  traced  directly  to  the  CHq  term  mentioned  above. 

The  parameters  CLa  and  Cp ^  were  varied  independently  of  C„a  in  order  to 

divorce  the  effect.  Variation  in  the  parameter  C.  were  investigated  se- 

a 

parately  and  in  combination  with  variations  in  CDa.  Increases  in  CLa 


alone  produced  mild  decreases  in  the  phugoid  frequencies  and  accompanying 
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increases  in  damping.  In  itself  this  is  not  particularly  important  but  it 
is  a  peculiarity  of  the  present  configuration  that  C0a  is  exactly  zero; 
thus,  it  w?s  feit  that  a  combined  variance  of  the  two  parameters  was  re¬ 
quired.  CD  variations  had  no  effect  on  the  characteristic  response  for 
a  large  number  of  CLa  values. 

The  parameter  was  assumed  to  be  negligible  in  the  preliminary 
analyses  because  of  the  high  pitch  rate  required  to  make  it  effective. 
Further  investigation  revealed,  however,  that  a  damping  term  was  present 
due  to  a  cable  whiplash  effect  following  a  displacement  of  the  vehicle. 

The  damping  occurs  because  the  energy  required  to  start  the  cable  oscil¬ 
lation,  which  of  course  is  damped  since  it  moves  upstream  along  the  cable, 
is  removed  from  tbe'vehicle  oscillation.  Therefore,  a  term  Cm^Q  was  in¬ 
cluded  in  the  computer  analyses  and  very  effectively  damps  the  short  period 
oscillation.  The  effect  of  Cm^g  on  the  phugoid  is  best  described  by  re¬ 
ference  fo  the  stabiiitv  boundary  of  Figure  18.  From  the  figure,  it  is 
seen  that  the  minimum  aliowatle  valde  of  Cma  for  stability  decreases  as 
Cfr^g  increases.  Th ■  s  effect  is  quite  similar  to  the  effect  of  C,,^  in  noi — 
mal  aircraft. 


The  effect  of  C-,^  variations  on  the  frequency  of  the  phugoid 
motion  was  significanr.  The  detailed  derivation  of  this  parameter  is  in¬ 
cluded  in  the  appendix;  only  note  for  this  discussion  that  the  large  values 
obtained  are  due  to  the  particular  configuration.  is  essentially  the 

restoring  force  of  the  svstem  but  also  appears  as  a  contributing  factor 
in  the  damping  of  the  phugoid  mode.  Increasing  ClT^  increases  the  damping 
but  decreases  the  frequency;  the  opposite  is  true  for  decreases  in  Cma- 
Along  with  these  changes  a  rather  significant  change  in  altitude  variation 
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for  a  given  input  was  noted,  an  increase  of  increasing  the  amplitude 
of  the  altitude  variation.  The  vehicle  can  not  be  compared  to  a  simple 
spring  mass  damper  case  for  illustrative  purposes;  the  system  is  simply 
too  complex.  contributes  to  both  the  restoring  force  and  damping  terms 

with  the  latter  effect  being  the  more  significant. 

The  cable  parameter  Cmc  is  not  actually  a  stability  derivative 
in  the  normal  sense  of  the  word.  Cmc  and  are  directly  related;  as 
one  increases  so  does  the  other  and  visa  versa.  This  may  be  seen  simply 
by  considering  the  equation  of  static  equilibrium.  As  discussed  pre¬ 
viously,  however,  even  though  is  not  a  stability  derivative,  it  plays 
a  major  role  in  determining  the  dynamic  response  of  the  vehicle,  thus  the 
inclusion  of  a  stabi  lity  boundary  depending  upon  as  Figure  19. 
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Y  LATERAL  DYNAMICS 

The  study  of  the  lateral-directional  dynamics  proceeded  in  the 
same  way  as  the  longitudinal  study  except  that  exclusive  use  was  make  of 
the  analog  computer.  The  study  was  again  specialized  to  the  basic  confi¬ 
guration  of  Figure  13  and  the  same  basic  flight  condition  of  170  mph  at 
sea  level  ar.d  a  cable  length  of  250  feet  was  uti  lized.  The  dynamic  sta¬ 
bility  derivatives  were  estimated  for  the  basic  flight  condition.  The  es¬ 
timation  procedures  are  described  in  Appendix  A. 

The  results  are  again  presented  in  the  form  of  stability  boundai — 
leu  of  the  important  parameters. 

Three  Dearee  of  Freedom  Studies 

These  studies  are  restricted  to  the  lateral-di rectiona I  modes  of 
motion.  Therefore,  the  terms  in  Equation  I  I  I -20  which  are  functions  of 
z  and  the  Nw  term  of  Equations  ill- 1 5  are  omitted.  In  other  words,  the 
tension  in  the  cable  is  assumed  to  remain  constant  throughout  the  lateral 
oscillations  and  the  body  is  assumed  to  be  perfectly  symmetrical  about 
its  centerline.  Then  substituting  the  approximation 

v  RS  U0P 

and  Equations  I  I  1-20  and  1 7 1 -22  into  Equations  I  1 1-1 5  after  replacing  the 
dimensional  stability  derivatives  with  the  non-dimensional  derivatives 
of  Table  I,  the  I atera I -di rect ional  eauations  of  motion  become: 


[D  -  (pUS/2n)Cyg  -  (To/mUbUje  -  [  (pSb/4m)CypD  +  (g/Uo) 

-  (T0  /  mUo )  ln  ]<p  +  [(dVJqHI  -  (pSb/4m)Cyr)  +  (To/mUoDly*  =  0  I-la 
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-  [(p^Sb/aixJC^  +  <a/lx)£TT0]e  +  [D2  -  <pUSb2/4lx)CifD  +  <a/ix>  • 

T0£r.]<P  ~  C (  ( lxz/ 1 +  <pUSb3/4lK)CAr02MI/Uo)  -  (a/ lx)  <T0/t)  ]y* 

=  0  3C-lb 

-  [(pU2Sb/2lz)Cnp  +  (e/lz)£7T0]e  -  [UX2/!X)D2  +  (pUSb2/4l2)CnpO 

-  (e/lz)T0jen]cp  +  [D3  -  (pUSb2/4lz)CnrD2)(l/Uo)  +  a0/\2)U7  -  l)}y* 

=0  5-1 c 


where  jJt  =  I  +  e/l ;  £n  =  I  +  a/n. 

Equations  3T-I  were  solved  on  the  analog  computer  to  obtain  the 
results  that  follow. 

The  lateral  equations  were  investigated  in  the  same  manner  as  the 
longitudinal;  first  variations  of  individual  parameters  then  combined  varia 
tions.  In  this  mode  there  are  long  period  and  short  period  osci I lations. 
The  dutch  roll  oscillation,  representing  the  short  period  mode  in  the  la¬ 
teral  case,  was  never  found  to  be  unstable  and  in  fact  little  could  be 
done  to  effect  the  damping.  The  mode  of  real  concern  is  not  the  dutch  roll 
however,  but  rather  the  pendulum  mode  becuase  it  is  the  mode  which  couples 
with  its  counterpart,  the  longitudinal  phugoid,  to  produce  an  unstable 
oscillatory  motion.  Thus,  the  following  remarks  concerning  parameter  vai — 
iations  wi il  be  oriented  towards  their  effects  on  the  pendulum  mode. 

Individual  variations  of  b'ng,  Cnp,  C^,  and  Cyr  produce  no 
appreciable  effect.  Decreasing  either  TD  or  C0r  decreases  the  damping  and 
frequency  of  the  pendulum  mode  while  increasing  either  has  the  opposite 
effects.  The  combined  variation  of  the  parameters  Cng  and  T0  has  a  neg¬ 
ligible  effect.  However,  decreases  in  both  Cnp  and  Cyg 


do  result  in  a 
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measurable  decrease  in  damping  and  frequency.  But  no  instabilities  were 
produced  by  magnitude  variations  of  these  parameters  within  the  limits 
of  fhe  computer. 

Two  combinations  did  produce  instabilities.  These  were  varia¬ 
tions  of  CnQ  with  C_gQ  and  Cr  with  Cn  .  Stabi  I  i ty  boundaries  were  deter- 
P  P  P  ' 

mined  for  these  variations  and  are  presented  in  Figures  20  and  21.  From 
Figure  20  it  is  seen  that  the  boundary  is  quite  similar  to  one  that  might 
be  obtained  for  a  normal  airplane.  The  oscillation  represented  is,  of 
course,  the  pendulum  oscillation  and  not  the  dutch  roll  for  which  no  in¬ 
stabilities  were  found.  Figure  21  indicates  that  below  a  certain  amount 
of  yaw  damping,  there  are  restrictions  on  the  minimum  value  of  Cng  for  a 
stable  system.  All  of  these  parameters  are  dependent  on  the  size  of  the 
vertical  tail  and  stability  will  in  general  be  increased  by  an  increase 
in  vertical  tail  size.  However,  the  location  of  the  added  area  with  re¬ 
spect  to  the  roll  axis  must  be  carefully  considered  so  that  the  value  of 
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21  LONGITUDINAL-LATERAL  COUPLED  DYNAMICS 

Following  the  investigation  of  the  longitudinal  and  lateral 
•nodes  separately,  the  equations  were  coupied  by  adding  the  z  dependent  and 
Nw  terms  of  Equations  1  I  I -70  and  I  I  I  -  I  5  to  Equations  Y- I  and  solving  the 
resulting  equations  simultaneously  with  the  longitudinal  Equations  TV-79. 

II  1-19  and  I  11-71. 

The  same  procedure  of  single  variation,  combined  variation  and 
finally  coupled  variations  of  the  significant  parameters  was  followed. in 
this  investigation.  Because  of  the  nature  of  the  coupling,  variations 
of  lateral  parameters  produced  no  changes  in  longitudinal  dynamics.  The 
opposite  was  not  true,  however.  For  example,  decreasing  Cm|_|Q,  which  de¬ 
creased  the  longitudinal  short  period  damping,  had  a  noticeable  effect 
on  the  dutch  roll  oscillation  of  the  lateral  mode.  This  i s  to  be  expected 
because  the  excitation  of  the  dutch  roll  is  from  the  yaw  moment  produced 
by  the  a  variation. 

Probably  the  most  significant  fact  to  come  to  light  in  this  in¬ 
vestigation  is  the  ex i stance  of  a  stray  coupling  between  the  phugoid  and 
pendulum  r..odes.  Aspects  of  this  coupling  are  discussed  in  the  following 
sections. 

Effect  of  Qr\a  and  Variable  Cable  Lengths 

Observations  made  of  the  vehicle  in  flight  indicated  that  the 
period  of  the  phugoid  was  an  inverse  function  of  the  cable  length.  For 
very  short  cable  lengths,  i .e.  five  feet  or  less,  the  vehicle  developed 
an  unstable  oscillation  coupling  the  longitudinal  and  lateral  modes  of 
motion.  The  analog  investigation  revealed  that  the  coupled  oscillations 
were  those  of  the  longitudinal  phugoid  and  lateral  pendulum  modes. 
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The  necessary  coupling  tern  was  found  from  wind  tunnel  data  to 
be  a  yawing  moment  due  to  change  in  angle  of  attack.  This  seems  to  arise 
from  the  use  of  a  drag  skirt  to  obtain  high  drag  coefficients.  The  result¬ 
ing  region  of  separated  flow  is  in  all  probability  sensitive  to  angle  of 
attack  changes,  thus  providing  the  condition  of  asymmetry  producing  Cna. 

As  was  discussed  in  the  section  on  Cable  Dynamics,  analytical 
studies  indicated  that  as  the  cable  became  shorter,  the  phugoid  frequency 
increased.  An  investigation  of  the  coupled  equations  of  motion  on  the 
analog  computer  bore  out  this  prediction  in  full.  Further,  as  the  plots 
of  cp  and  z  (Figure  22),  which  represent  the  pendulum  and  phugoid  modes 
respectively,  indicate  the  coupled  oscillations  do  indeed  develop  and  in¬ 
crease  in  intensity  as  the  cable  becomes  shorter. 

Increases  in  phugoid  frequency  aggravate  the  coupled  motion  dis¬ 
cussed  above.  The  most  significant  parameter  in  this  respect  is,  of  course, 
rope  length,  although  others  do  exist.  Decreasing  and  the  combination 
of  decreasing  Cn^  and  simultaneously  increasing  C,^  have  relatively  strong 
influences  with  respect  to  increasing  the  phugoid  frequency.  In  the  la¬ 
teral  mode,  any  increase  in  frequency  of  the  pendulum  mode  was  accom¬ 
panied  by  an  increase  in  damping  with  the  net  result  being  neglibible. 

Besides  configuration  change's  of  the  vehicle,  cable  character¬ 
istics  were  varied  to  determine  if  they  might  have  any  influence  on  this 
aspect  of  the  study.  As  can  be  seen  from  Figure  23,  a  plot  of  phugoid 
period  vs.  length,  the  effect  of  varying  a°,  essentially  the  weight  to 
cross  section  ratio  of  the  cable,  decreases  rapidly  with  the  smaller  values 
of  l.  Thus,  it  is  felt  that  a°  is  relatively  unimportant. 
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APPENDIX  A 

Derivation  of  Stability  Derivatives 

The  longitudinal  stability  derivatives  were  obtained  from  the 
first  wind  tunnel  report  from  DTM3;  however,  several  of  the  lateral  sta- 
biliry  derivatives  were  calculated  analytically,  a  discussion  of  the 
technique  follows. 


Let  the  distance  from  the  c.g.  to  the  aerodynamic  center  of  the 
vertical  stabilizer  be  e.  The  yaw  velocity,  r,  produces  a  side  velocity 
at  the  tai  I . 


VT  =  er 

The  si de  slip  angle  at  the  tail  is  then, 

3r  =  -  VT/lio  =  -  er/lio  ~  radians 

The  lift  curve  slope  of  the  vertical  tail  is  taken  as  CL  ,  then  defining 

Pt 

Cvo  —  Ci  .  ST/S,  the  total  force  due  to  sideslip  becomes: 

“t  Pt 

Ft  =  (CygT  PT )  ( 1/2  pUsS) 

The  yawing  moment  due  to  this  force  is, 

N  =  -  eFT 


Therefore, 


C  / 

CN  =  +  e  YPt  •  er/ UU  -  1/2  PU3S 

1/2  pU3Sb 


and  after  cancellation  and  multiplication  by  b/b 


CN  =  (2es/bs)Cv 


rb/2U 


'N  -  ^ypT 

Differentiating  this  expression  with  respect  to  rb/2U,  yields 


yPr 


.  Cnr  as: 


Cnr  = 

3rb 

2U 


2(e/b)2C 
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The  coefficient  Cyr  may  be  obtained  in  much  the  same  manner. 
From  the  equation, 

cy  =Cy^(-er/U0)  =  -  (2e/b)  Cyp  (rb/2U0) 

Cyr  is  obtained  by  differentiating  this  expression  for  Cy  with  respect 
to  rb/2Uo,  thus. 


C 


V- 


Cyr  =  -  (2e/b)CygT 


(See  Figure  Al) 
Defining  vp  =  f 

a  =  w/U 


=  ao 


+  Aa 


UL  =  U0  +  Ur 
UR  =  Ub  -  Ur 
oic  =  w/U0 

therefore,  w/U  -  w/U0  =  Aa 

The  incremental  changes  in  angle  of  attack  on  the  left  and  right 
wing  are  respectively. 


AaL  =  -  do  (Ur/Uo-HJ,0 


and 


AaR  =  +  ao  (Ur/Uo~Ur) 

From  the  foregoing  analysis,  AL„  takes  the  following  form 

ALr  =  Lr  -  Uo  =  1/2  pUr2  S  CLn  ar  -  1/2  pUo3  S  CL  ao 

w  a 


which  becomes 

ALr  =  1/2  PS  CL  [(Uo  ~  Ur)2(Oo  +  Aar)  -  (Uo2Oo)] 
This  is  further  reduced  to  the  final  form  of 


ALr  =  1/2  pS  Ct  (-  UpUf-do ) 
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The  same  procedure  was  followed  in  obtaining  ALl  , 

ALl  =  1/2  pS  CL^  UoUrOo 

Defining  AL  as  being  positive  in  the  negative  z  direction,  the  expression 
for  the  ro I i  moment  i  is: 

i  =  -  gALr  +  gALt 

which  after  substitution  and  cancellation  reduces  to, 

&  =  "f  3pSC|_^  UoUrOo. 

where  g  is  the  distance  from  the  centerline  to  the  M.AvC. 

Thus,  becomes, 

Ctj  =  (4g/b2)  CLa  do  (Urb/2Uo ) 

si  nee 


Ur  =  vr  cos  0  =  rf  cos  0 

and  f  =  f  cos  Q 

Upon  substitution,  Ur  becomes: 


Ur  =  gr 

Thus,  in  final  form,  C i  is 

CA  =  (2g/b)a  Cla  ao  (rb/2Uo) 

By  differentiating  with  respect  to  rb/2U0  and  making  the  substi¬ 
tution  of  do*  =  do  "do  ,  C a  is, 

i  *r 

CAp  =  (  9  )3  C  do* 
b/2 


Cyp: 

Due  to  the  symmetry  of  the  vertical  stabiilzers  about  the  x  axis, 
the  contribution  of  these  surfaces  is  negligible.  Also,  due  to  the  small 
perturbation  assumption,  the  contribution  of  the  horizontal  surfaces  may 
be  taken  as  zero.  Therefore, 

Cyp  =0 
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% 

% 

at  the  m.a.c. 
Agai n. 


has  been  estimated  by 
resulting  from  the  ro 


calculating  the  average  angle  of 
I  velocity  p.  (See  Figure  2A.) 


attack 


t  =  -  gALr  +  gALL 


where 


ALr  =  CL  ar  gS;  ALl  =  CLa  aL  qS 

and 


ar  sa  +  gp/D0;  aL  -  gp/Uo 
Thus,  upon  substitution,  l  becomes, 

l  =  -  2g  CLa  gS  gp/UQ 
wh  i  ch,  when  nondimensional  ize<4  is: 


C£  =  -  (  9  )a  CL 

b/2 

By  differentiating  with  respect  to  pb/2U0, 

C,  =  -  (  _9_  ) 
P  b/2 


pb/2Uo 

i s  found  to  have  the  form, 

2  r 


The  vertical  fins  again  contribute  nothing  to  this  parameter  be¬ 
cause  of  their  symmetrical  location.  The  wing  contribution  stems  from  two 
sources;  the  change  in  drag  due  to  the  new  angle  of  attack  resulting  from 
roll  and  the  moment  produced  as  a  result  of  the  non-parallel  alignment  of 
the  I i ft  vectors. 

Because  only  small  perturbations  from  the  equilibrium  position 
where  C0a  is  zero  are  considered,  the  contribution  from  the  first  source 
i s  neg I i gi b le. 

Establishing  the  positive  value  of  AD  as  pointing  in  the  nega¬ 
tive  x  direction  the  expression  for  the  yaw  moment  N,  becomes,  (See  Figure 


A3) 
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N  =  gADr  -  gADL 

From  Figure  4A,  it  may  be  seen  that 

ADr  =  -  Lr  sin  Aar  -  LrAar 


and 


ADl  =  -  Ll  sin  AaL  =  -  LLAaL 


The  expression  tor  Ll  i s  as  follows. 


a  I  sc 


Ll  =  1/2  pSUl2CLl  =  1/2  pSUL2CLa  aL 


Lr  =  1/2  pSUr2CLa  ar 
6y  making  the  identification, 

a  =  Oq  +  Aa,  =  (wq/Uq)  +  (w/U) 

and  substi tjt : ng  the  expression  for  a  into  the  previously  developed  terms 
for  Lr  and  Ll ,  they  become, 

Lr  =  1/2  pS‘Jr2CLa[(w0/U0)  +(wr/Ur)  ] 
a.'.d 

Ll  =  1/2  pSUr2CLa[(w0/LI0)  +  (WL/UL)] 

Then  with  wL  =  -  gp  and  wr  =  +  gp,  the  expressions  for  ADr  and 
AOl  become, 

ADr  =  -  1/2  pSUraCla[oo  +  (gp/Ur) ]Cgp/Ur) 

and 

ADl  =  -  1/2  PSUL2CL  [cto  -  Cgp/UL)]<-  gp/'Jr) 

Using  the  identities, 

UL2  =  Up2  +  wL2 ;  Ur2  =  U02  +  wr2 

and  the  inequalities, 

w r2  «  Uc3;  wLs  «  U02 
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the  expression  for  N  becomes, 

N  =  1/2  pSg  [2U0CLa  Oq  gp] 

Then, 

Cn  =  -(4g2/b3)  CLa  cto  (pb/2Uo) 


Differentiating  with  respect  to  pb/2Uo ,  the  desired  derivative 
i  s  obtai ned  as: 


(b/2) 


3  °La 


The  values  of  the  longitudinal  and  lateral  stability  derivatives 
for  the  basic  configuration  at  a  velocity  of  170  mph  at  sea  leval  are  pre¬ 
sented  i n  Tab le  A- I . 


TABLE  A- I 

NOND I  MENS  1 ONAL  STABILITY  DERIVATIVES 
BASIC  CONFIGURATION 


ye 

-  1.605 

C*r 

=  + 

0.02 

ne 

= 

+  0.647 

CyP 

=  0 

H 

= 

-  0.067 

°nP 

=  - 

0.02 

Yr 

= 

+  2.86 

C  i 
£P 

=  - 

0.83 

nr 

= 

-  2.37 

=  - 

6.88 

mq 

= 

-  0.076 

CL<x 

= 

3.58 
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APPENDIX  B 

Calculation  of  Equi librium  FI i Qht  Conditions 

The  equilibrium  equations,  for  tho  case  of  zero  flight  path  an¬ 
gle,  are  obtained  from  Equations  I  11-13  as: 

Xo  +  TXo  =  0 

Z0  +  W  +  TZq  =  0  B-l 

+  T*p  =  0 

The  cable  tension  components  in  these  equations  are: 

Txo  =  To  cos  X0 

=  -T0 si n  \q  B-2 

THo  =  -Tc[e  sin  (oo-X0)  +  a  cos  (oo-Xq)] 

The  aerody.nami  c  terms  are  then 

X0  =  -  1/2  PU2S[CDmin  (CLas  (ao  -OoL  )/nA)  ] 

Z0  =  -  1/2  pU2S  CUa(ao-aoL)  B-3 

H>  *  1/2  PU8Sc[CM0l  +  C„a  Cao-tXoL)] 

where  the  subscript  0l  indicates  the  zero  lift  condition.  Substitution 
of  Equations  B-2  and  B-3  into  Equation  B-l  and  assuming  ao  «  XD ,  yields 
the  equations: 

T0  cos  X0  -  1/2  plt’sCCo  Cigs  (gQ-aoL)3  ]  =  0  B-4 

TTA 

T0  sin  Xo  -  W  +  1/2  PU3S  CLa  (ao'  OoL)  =  0  B-5 

1/2  pUsSc  (CHQ  -  CH  ao  )  +  e  T0  sin  Xo  -  aT0  cos  X0  "  [T0(a  cos  X0 

l  si  L 

+  b  sir.  X0)  -  1/2  pU2Sc  CMa]ao  =  0  B-6 


These  three  equations  must  now  be  solved  for  the  unknown  factors  T0,  X0 
and  cxq. 
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Substi tut  ion  of  Equations  B-4  and  B-5  into  B-6  then  gives  the 
cubic  equation  for  Oq 

Oq3  -  (2ool  +  (anA/eC^)  -  a/e)ao8  +  [ (Cp^ j  ^ttA/Cl^2)  t-Oo^2 
+  (WnA/l/2  pUaSCLa2)  a/e  +  (3nA/eCLa)  -  c/e  )  +  nA/CUa 

-  2  a/e  do  ]ao  *  (WttA/I/2  pliaSCLa2)  -  (ttA/C,^)^  +  a/e  (Cqm(  NnA/CLa2) 
+  a/e  0oL2  -  (CnA/eC^HC^  -  C^Oo^  =  0  B-7 

Equation  B-7  is  ihen  solved  numerically  for  Oq . 

With  ao  known.  Equations  B-4  and  B-5  are  then  solvea  for  \0 

X0  =  tan*1  W  +  _  B-8 

q5[CQj^|jg  +  CCl^2/ttA)  (aoL8+  do2-  2aoLao) 

The  cable  tension  is  then  obtained  from  Equation  B-5, 

T0  =  W  ~  clSCLq^ao  ~  aol  >  B-° 

si  n  \0 

TABLE  B-l 


u 

<k> 

x0 

To 

alt. 

150 

-0.94 

60.25 

70.0 

$L 

170 

-1.35 

41.65 

95.5 

SL 

200 

-1.96 

33.95 

1  12.4 

SL 

250 

-2.45 

26.60 

181.5 

SL 

300 

-2.77 

21.90 

252.0 

SL 

350 

-2.94 

19.35 

345.0 

SL 

170 

-1.29 

43.20 

93.0 

2000' 

170 

-1  .  10 

53.00 

78. 1 

5000' 

250 

-2.42 

27.55 

172.7 

2000' 
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APPENDIX  C 

Computer  Analysis 

Early  if  the  investigation,  it  became  apparent  that  the  problem 
would  lend  itself  nicely  to  analysis  by  analog  computer  techniques.  The 
computer  work  progressed  concurrently  with  analytical  studies  following 
the  standard  approach  of  first  examining  the  uncoupled  equations  of  motion 
and  then  a  complete,  i.e.  coupled,  set  of  equations. 

The  longitudinal  equations  are  really  only  three  in  number,  the 
highest  order  being  second.  What  may  appear  as  a  fourth  equation  is  the 
integral  expression  for  z  which,  when  multiplied  by  the  proper  coefficients, 
becomes  the  cable  contribution  to  the  system  dynamics.  Electronic  multi¬ 
pliers  were  used  since  they  are  standard  equipment  on  the  Pace  TR-IO  com¬ 
puters  used  for  this  work. 

The  lateral  equations  were  three  in  number,  the  highest  order 
bei ng  thi rd. 

The  coupled  equations  were  six  in  number,  and  are  essentially 
the  same  equ.rMons  as  were  treated  independently,  the  only  significant 
change  being  the  addition  of  the  Cno  term  in  the  lateral  equation.  This 
is  discussed  in  another  section.  The  coupling  required  the  addition  of 
three  more  multipliers  and  a  number  of  amplifiers  and  pots.  When  com¬ 
pletely  coupled,  three  units  were  slaved  employing  some  52  operational 
amplifiers,  a  tribute  to  the  flexibility  and  capability  of  the  Pace 
equi pment. 

Figure  C-l  is  a  schematic  of  the  wiring  diagram  used  in  the 


computer  study. 


FIGURE  L  COMPLETE  SYSTEM  REPRESENTATION 


FIGURE  2.  DISTURBED  CABLE  ELEMENT  IN  VACUUM 


FIGURE  3.  CABLE  ELEMENT  IN  EQUILIBRIUM 


FIGURE  4.  CONFIGURATION  OF  CABLE 

SUPPORTING  A  ZERO  DRAG  BODY 


FIGURE  &  FORCES  ON  A  CABLE  IN  A  UNIFORM  STREAM 


FIGURE  12.  GEOMETRY- DEPENDENCE  OF  X  VARIATION  ON 

LINE  LENGTH 


Wing  area  S  -  5.5  ft 
Wing  span  b  -  4  ft 
Mean  chord  c  =  133  ft. 

a  =  ,625  ft 
e  =  2.67  ft. 
g  =  ,964  ft 
Aspect  Ratio  A  =  2.91 
Root  chord  Cr  =  1.5  ft 
Tip  chord  Ct  *  1.25  ft 
Sweepback  angle  jv.  =  37 


FIGURE  13.  BASIC  BIRD  CONFIGURATION 


I 

I 


FIGURE  16.  APPROXIMATE  ZERO-POLE  LOCATIONS  FOR  THE 
TRANSFER  FUNCTIONS  OF  EQUATIONS  IE  -  14 


FIGURE  17.  GEOMETRY  FOR  STATIC  EFFECT  ON  VARIABLE  TENSION 


FIGURE  Cl.  ANALOG  DIAGRAM 


